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Introduction 



La physique des particules est apparue pour traiter rinfiniment petit et chercher 
a comprendre I'organisation de la matiere a son niveau le plus fondamentale. L'as- 
trophysique a son tour s'ouvre a I'infiniment grand et tente de percer les mysteres 
de I'organisation de I'univers dans son ensemble, de meme, la cosmologie s'interesse 
a rinfiniment ancien en essayant de reconstituer I'histoire de I'univers depuis ses 
premiers instants. 

Autour de ces trois filieres s'est developpee une nouvelle discipline, I'astroparticule. 
Un domaine vaste et fortement interdisciplinaire, ayant pour objectif principale 
I'etude de toutes particules ou rayonnement venant de I'espace, la mise en evidence 
et la comprehension des phenomenes cosmiques de haute energie. II s'agit d'etudier 
(ou decouvrir) des objets au sein desquels ont lieu des processus non thermiques, 
capables d'emettre des particules (notamment des rayons cosmiques (RC)) de haute 
energie au dessus du GeV. 

En effet la radiation cosmique, mis en evidence par Victor Hess en 1912, est une 
composante principale de la galaxie, vu que sa densite d'energie est comparable, si 
elle n'est pas superieure, a celle des autres radiations presentes dans I'univers. 
Les experiences souterraines jouent un role important dans I'etude des differentes 
caracteristiques des RC telles que la composition chimique, les mecanismes d'accele- 
rations, I'astronomie des rayons gamma et des neutrinos en plus de la recherche des 
phenomenes rares (desintegration du proton, monopoles magnetiques, nuclearites, 
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etc.). 

L'experience MACRO, situee en Italie, constituait I'une des plus grandes expe- 
riences souterraines. Elle a permis de couvrir un vaste domaine de recherche dont 
on cite : la recherche des particules rares tels que les monopoles magnetiques |38| et 
les nuclearites [86], I'etude des neutrinos atmospheriques [T][2] et de leur oscillations 
[3], I'etude des muons de haute energie du bas [4] [55] et I'etude de la composante 
penetrante des RC. 

Parmi les questions de recherche qui presentent un interet fondamental, on dis- 
tingue I'etude de I'anisotropie temporelle du flux de muons cosmiques. Des mo- 
dulations dans les distributions temporelles du temps d'arrivee de muons out ete 
observees par quelques experiences. Deux d'entre elles out mis en evidence la pre- 
sence d'un signal module avec une periode de 4.8 h, exactement egale a la periode 
orbitale de Cyg X-3 dans la region des X du spectre. D'autres ont observe des mo- 
dulations dans les distributions du temps d'arrivee des muons dues a la presence 
d'une composante qui n'est pas aleatoire. Les resultats de ces experiences n'ont pas 
ete conflrmes par d'autres groupes et le probleme reste confus. 

Dans ce travail on se propose d'etudier les variations du flux des muons cos- 
miques, et ceci en utilisant les donnees collectees par le detecteur souterrain MA- 
CRO. Nous disposons d'un grand nombre de muons (environ 38 millions d'evene- 
ments) detectes durant une periode de plus de 6 ans. 

Cette these est organisee en 5 chapitres. Dans le premier, nous passons sur un 
rappel sur les rayons cosmiques, les sources et les mecanismes d'accelerations, leurs 
interactions dans I'atmosphere ainsi que la perte d'energie des muons dans la roche. 
Dans le second, nous allons donner une vue sur quelques experiences souterraines 
et notamment l'experience MACRO. 

Le troisieme chapitre sera consacre a I'etude des distributions des temps d'arrivee 
des muons cosmiques avec une energie plus grande que 1.3 TeV au sommet de 
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la montagne de Gran Sasso et collecte durant la periode allant de 1995 jusqu'a 
2000. Un grand echantillon de muons a ete utilise compare a celui utilise dans les 
travaux ulterieurs de MACRO. Cette etude nous permet de chercher I'existence de 
modulations dans le flux de muons cosmiques et de conclure sur la nature de la 
distribution des temps d'arrivee des RC afin de se situer par rapport aux resultats 
obtenus par les experiences. 

Dans le quatrieme chapitre, nous presentons les resultats de I'etude des varia- 
tions du flux de muons en utilisant deux approches complementaires : 

- Recherche de clusters d'evenements avec la methode Scan statistics. 

- Recherche de variations periodiques avec la methode de Lomb-Scargle. 

La derniere partie de ce travail concerne I'etude de la perte d'energie des nuclearites 
dans difl^erents milieux, notamment dans I'atmosphere, afln d'evaluer la contribution 
de I'atmosphere dans le processus de perte d'energie. 
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Chapitre 1 



Generalites sur les rayons cosmiques 



1.1 Introduction 

Le rayonnement cosmique demeiire un des problemes centraux de I'astrophy- 
sique, suscitant de nombreuses questions bien souvent interconnectees. Son etude 
permettra d'etablir des liens tres riches entre les differentes parties des sciences de 
I'univers, touchant notamment a I'etude du milieu interstellaire, des champs ma- 
gnetiques ou de I'interaction des particules avec les plasmas, a la modelisation des 
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sources energetiques telles que les noyaux actifs de galaxie ou les sursauts gamma, a 
I'astronomie gamma, la formation des etoiles, a la nucleosynthese spallative, etc. En 
effectuant plusieurs vols en ballon avec des electroscopes, le scientifique Victor Hess 
apporta la preuve, en 1912, que ce rayonnement vient de I'espace. Ce rayonnement 
fut surnomme "rayonnement cosmique" par Robert Millikan en 1925. C'est grace 
a cette decouverte que la physique des particules experimentale prit son essor et 
son etude a conduit entre 1930 et 1950 a la decouverte du positron, du muon, des 
mesons n et k ainsi que les baryons E, A et S. 

En 1938, Pierre Auger et ses collaborateurs observerent des coincidences en temps 
entre plusieurs detecteurs de particules separes de plusieurs metres. Auger emis 
alors I'hypothese selon laquelle les particules detectees au sol etaient les consti- 
tuants de gerbes atmospheriques initiees par des rayons cosmiques dont I'energie 
pent aller au dela de 10^^ eV. 

Aujourd'hui, les observations d'evenements autour de 10^° eV restent encore une 
enigme. Une partie de la comprehension des rayons cosmiques d'ultra haute energie 
passe par I'etude de leurs interactions dans I'espace interstellaire et I'atmosphere. 
Majoritairement constitue de particules chargees, le rayonnement cosmique frappe 
I'atmosphere terrestre de maniere continue. Le flux de particules chargees qui 
heurtent I'atmosphere pent atteindre 1000 particules/m^s. 
Les rayons cosmiques peuvent etre classifles en trois types. 

- Les RC primaires charges, sont constitues d'environ 90% de protons, 9% de 
particules a et le reste par des noyaux atomiques tels que (C, O, N, Mg, Fe,...). 
Toutefois, cette composition varie avec I'energie, par exemple, a des energies 
de I'ordre de 400 TeV on a environ 12% de protons, 25% de particules a, 
26% de CNO, 15% de Si-S et 21% de noyaux avec Z>17. 

- Les RC secondaires, sont constitues de mesons vr, K et des p,n, etc., re- 
sultant de I'interaction des primaires charges avec des noyaux atomiques de 

il GeV= lO^eV, 1 TeV= lO^^eV, 1 PeV= lO^^eV et 1 EeV= lO^^eV. 
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I'atmosphere terrestre a une altitude moyenne de I'ordre de 20 km. 
- Les RC tertiaires, sont constitues de muons, neutrinos, photons, gerbes elec- 
tromagnetiques et des electrons resultant de la desintegration des mesons 
charges et neutres des RC secondaires. 

Du point de vue energetique les rayons cosmiques representent une composante fon- 
damentale de I'univers. Leur densite energetique est approximativement 1 eVjcrr? 
comparee a la densite d'energie de la lumiere stellaire qui est de I'ordre de 0.6 
€\l jcrr?. On dispose alors, d'un veritable echantillon de la matiere extraterrestre. 
Si aux plus basses energies, la propagation des rayons cosmiques charges dans les 
champs magnetiques va uniformiser leurs directions d'arrivee, rendant ce rayon- 
nement isotrope autour de la terre, aux plus hautes energies, ceux-ci vont avoir 
une rigidite suffisante pour permettre de "pointer" leur source. Cette isotropic de 
la direction d'arrivee rend difficile la localisation des sources et par consequent il 
est difficile d'etablir I'origine galactique ou extragalactique des rayons cosmiques 
primaires. L'absence de toute anisotropic dans la direction de provenance pent etre 
due a une distribution isotrope des sources, aux effets de propagation dans I'espace 
interstellaire ou a un manque de statistique. En effet, en traversant I'espace entre la 
source et la terre, les RC peuvent interagir avec les photons du fond cosmologique 
et le gaz interstellaire et peuvent subir egalement des deviations provoquees par les 
champs magnetiques galactiques et extragalactiques. 

La problematique des rayons cosmiques se presente a la fois sous des aspects theo- 
riques, phenomenologiques et observationnels. Des travaux mettent I'accent sur 
I'acceleration des rayons cosmiques, tandis que d'autres se focalisent sur leur pro- 
pagation, ou sur leurs interactions avec le milieu environnant, depuis leurs sources 
jusqu'a I'atmosphere terrestre. 

La comprehension des phenomenes a I'origine de la production des rayons cos- 
miques, les mecanismes d'acceleration et de propagation dans le milieu interstel- 
laire et intergalactique, necessite la connaissance de leur composition chimique, du 
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spectre energetique, de la distribution du temps et de la direction d'arrivee. 
Pen de questions trouvent leurs reponses sur I'origine des RC, cependant on salt 
clairement que ceux de haute energie proviennent de I'exterieur du systeme solaire. 
On ne connait pas encore ni les sources responsables de remission des RC ni les 
mecanismes qui favorisent leur acceleration a de tres hautes energies (> 10^^ eV). 



1.2 Spectre energetique des RC 

Plus de 90 annees d'observation des rayons cosmiques ont permis de construire 
leur spectre en energie (figure (ll.ip ) qui s'etend de maniere remarquablement conti- 
nue sur environ 13 decades en energie (entre ~ lO^eV^ et ~ 10^^) et 32 decades en 
flux (entre quelques milliers de particules par m^ et par seconde aux basses energies 
et de I'ordre d'une particule par km^ et par siecle pour les plus energetiques) . 

Get intervalle etant tres large, il parait clairement que pour le couvrir experi- 
mentalement, I'utilisation de diverses techniques de detection est exigee, du petit 
detecteur tel que la chambre a emulsion dont I'acceptance est de I'ordre de 0.1 
rn? aux grands detecteurs dont la surface pent atteindre des km?. Tout cela est 
egalement accompagne de nombreux problemes que ce soit au niveau de I'analyse 
des donnees (normalisation des resultats produits par les difl'erents detecteurs) ou 
au niveau de la technique experimentale (detection des particules de tres hautes 
energies). 

En efl'et le flux difl'erentiel pent s'exprimer comme une simple loi de puissance [5] : 

dN ^ ^ , , 

-ocE- (1.1) 

ou E est I'energie par nucleon de la radiation, 7 est I'indice spectral qui depend de 
I'energie et N est I'intensite des RC. 

En efl'et, le spectre n'est pas un pur spectre de puissance, c'est a dire I'exposant 
7 depend de I'energie. Pourtant, dans une tres large bande d'energie il s'avere que 
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Fig. 1.1 Spectre des rayons cosmiques (flux en fonction de I'energie). 
En dessus sont indiquees approximativement les regions couvertes par 
les mesures directes et indirectes. 
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I'equation (ll.ip donne une bonne approximation du spectre. Comme on pent le voir 
sur la figure (ll.ip il existe plusieurs regimes de ce spectre a savoir : 

- En dessous de quelques MeV le fiux est domine par les particules solaires. 
Celles-ci constituent le vent solaire qui infiuence les rayons cosmiques d'energie 
allant jusqu'a quelques GeV par le biais de la modulation solaire. 

- Entre ~ 0.1 GeV et 100 TeV I'indice spectral est 7 ~ 2.7. Ceci est assez bien 
explique par les mecanismes d'acceleration classiques. II faut cependant tenir 
compte du fait que jusqu'a une energie de quelques GeV, la loi de puissance 
est modifiee par la modulation solaire. 

- Entre 10^"^ eV et 10^^ eV I'indice spectral devient 7 ~ 3.2. La zone situee 
autour de 3-5x10^^ ou a lieu le changement de la pente est appelee genou 
(Knee). Dans cette gamme d'energies les fiux sont tellement faibles (de I'ordre 
d'une particule par m^ et par an) que I'on ne pent plus faire les mesures de 
maniere directe. On doit alors observer les cascades atmospheriques. 

- Une autre rupture de pente apparait a 10^^ eV communement appelee la che- 
ville (ankle). Celle ci marque I'entree dans le domaine des rayons cosmiques 
d'ultra haute energie pour lesquels I'indice spectral se radoucit et devient 
7 ~ 2.8. Ici les fiux sont tellement faibles (de I'ordre de 1 particule par km^ et 
par siecle) que les resultats experimentaux sont rares. lis ne permettent pas a 
I'heure actuelle d'eclaircir le mystere de I'origine de telles particules. En effet 
il parait irrealiste de produire de telles energies dans les mecanismes d'ac- 
celeration classiques. Une solution pourrait venir d'objets extragalactiques 
comme les noyaux actifs de galaxies. Cependant, I'interaction de particules 
avec une telle energie avec les photons du fond cosmologique (CMB) devrait 
provoquer une perte d'energie par photoproduction de pions. C'est la cou- 
pure appelee GZK [6]. L'experience AGASA [8] semble voir une anisotropic 
et un exces d'evenements au dessus de la coupure GZK par rapport a une 
source extragalactique, ce qui n'est pas le cas du detecteur HiRes [9j. En 
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tenant compte des erreurs systematiques importantes, leurs mesures restent 
cependant compatibles et ne fournissent done pas d'indice significatif. Pour 
commencer a lever la voile il faudra attendre d'avoir une bonne statistique 
avec I'experience AUGER [TO] . 
Ces changements peuvent etre dus a des modifications substantielles dans les me- 
canismes de production et d'acceleration des RC de haute energie, on etre lies a 
la propagation du rayonnement cosmique dans I'espace. Cette derniere idee a ete 
developpee par le modele de diff'usion des RC dans la galaxie dit "leaky box" [TI] . 
Selon ce modele, la galaxie est consideree comme un reservoir confinant des RC 
qui permet a ces derniers d'avoir des probabilites non nulles de s'enfuir. En eff'et, 
les particules produites dans notre galaxie traversent le milieu interstellaire avant 
d'atteindre la terre. Ce milieu est constitue de nuages de gaz neutre et ionise et 
il est caracterise par un champ magnetique galactique desordonne B (~ 3 x 10~^ 
gauss). Une particule de charge Z qui traverse ces nuages magnetiques subit des 
defiections et aura un rayon de courbure, 

R[cm] = E[eV]/300 ■ B[gauss]Z (1.2) 

Lorsque I'energie augmente, le rayon augmente, ce qui signifie que beaucoup de pro- 
tons (Z=l) atteignent des rayons superieurs an rayon de la galaxie elle meme. Cet 
eff'et tend a les faire disparaitre dans I'espace intergalactique et ils ne peuvent plus 
par consequent contribuer au ffux de RC. Naturellement il y a d'autres particules 
qui arrivent d'autres galaxies. 

A leur tour, les elements chimiques simples qui constituent les RC suivent Fai- 
lure de la fonction exponentielle. Sur la figure ( II. 2p [l2j . on presente le spectre 
energetique individuel des noyaux presents dans les RC et qui pent etre decrit ega- 
lement par une loi exponentielle dNi/dE = KE~'^^ avec 7j et K qui changent selon 
que I'energie soit superieure ou non a I'energie du genou. 
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Fig. 1.2 Spectre energetique individuel de I'hydrogene, Vhelium, du car- 
hone et du fer. 
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1.3 Composition chimique des RC primaires 



L'etude de I'abondance relative des elements du systeme solaire (basee sur I'ana- 
lyse des materiaux terrestres et des meteorites et les mesures spectroscopiques des 
etoiles tel que le soleil) et celle des isotopes presents dans les rayons cosmiques per- 
met d'obtenir des informations importantes sur la nature des sources emettrices et 
sur le mode de propagation de ces particules dans le milieu interstellaire. L'abon- 
dance chimique des elements dans les RC est connue au moins jusqu'a quelques 
dizaines de TeV/nucleon a partir des observations perforees au niveau de la haute 
atmosphere a I'aide des ballons, des satellites,...). On pent voir sur la figure (II. 3p 
les abondances relatives des diff'erentes especes dans le rayonnement cosmique com- 
parees a celles du systeme solaire [12] ; les cercles pleins sont relatifs aux donnees 
a basse energie (70 - 280) MeV/A, les cercles vides sont relatifs aux donnees des 
rayons cosmiques a haute energie (1000 - 2000) MeV/A et les losanges pour I'abon- 
dance isotopique dans le systeme solaire |13j. 

On observe une similitude entre les deux compositions chimiques a quelques ano- 
malies notoires pour certains noyaux. Des differences apparaissent cependant pour 
les elements immediatement plus legers que le fer. Le meme phenomene touche 
egalement les elements (Li, Be et B) pour lesquels les abondances cosmiques sont 
beaucoup plus importantes. Ce phenomene pent etre interprete comme I'existence 
de reactions de spallation dans le milieu interstellaire. En effet, la fragmentation 
d'elements comme le carbone, I'azote ou I'oxygene (CNO) sur des protons du milieu 
interstellaire produit une forte quantite d'elements plus legers comme le Lithium, 
le Beryllium et le Bore enrichissant leurs portions dans le rayonnement cosmique. 
La comparaison des abondances entre le CNO et le Lithium, le Beryllium et le Bore 
nous renseigne sur la quantite de la matiere interstellaire traversees. 
II existe aussi un deficit de I'hydrogene et de I'Helium dans le rayonnement cos- 
mique par rapport a I'abondance stellaire. Ceci pent etre du a la difficulte d'ioniser 
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Fig. 1.3 Abondance relative des elements dans le rayonnement cosmique 
(ligne solide) et dans le systeme solaire (ligne en pointillee) normalise 
a 100 pour le Si. 
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ces elements, les rendant moins disponibles pour I'acceleration. On constate aussi 
la presence de I'effet d'appariement dans les deux compositions. En effet les noyaux 
pair-pair sont en general stables et plus abondants. 

Certaines de ces differences peuvent trouver des explications dans les modeles de 
propagations des RC dans le milieu interstellaire. Avant d'atteindre la terre, cer- 
tains noyaux subissent la fragmentation et produisent des noyaux plus legers. Ce 
qu' on observe localement est constitue d'une composante primaire provenant direc- 
tement de la source et d'une composante secondaire resultant de la fragmentation 
de certains noyaux apres interaction avec d'autres particules durant le passage dans 
le milieu interstellaire. De ce fait, les elements (Li, Be, B) et (Sc, Ti, V, Cr, Mn) 
peuvent resulter de la fragmentation du (C, O) et Fe respectivement. Cette ca- 
racteristique temoigne du role important de la transformation de la composition 
chimique des RC lors de leur propagation a travers I'espace interstellaire et peut- 
etre au sein des sources (dans la region de la production ou de I'acceleration des 
rayons cosmiques). 



1.4 Direction d'arrivee 

A faible energie, les rayons cosmiques subissent des deflexions de leurs trajec- 
toires causees par les champs magnetiques galactiques et extragalactiques. 
Au dessus de lO^^eV les rayons cosmiques sont peu devies de leur trajectoire par 
les champs magnetiques galactiques et extragalactiques (de I'ordre de quelques de- 
gres), leur direction d'arrivee devrait done pointer vers la source. 
Autour de lO^^eV, les experiences AGASA et Fly's Eye ont reporte une anisotropic 
en provenance du centre galactique avec un deficit correspondant dans la direc- 
tion oppose [IS] [21]. Au dela de lO^^eV, les resultats des differentes experiences 
semblent s'accorder sur une repartition isotrope des rayons cosmiques. AGASA a 
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notamment detecte une telle repartition avec la presence de trois doublets et d'un 
triplets d'evenements dans des cercles de diametres 2.5°. Aucun objet astrophysique 
n'a ete correle avec la direction de ces multiplets (figure (ll.4l) )[22]. 




Fig. 1.4 Directions d'arrivee, en coordonnees equatoriales, des rayons 
cosmiques au-dessus de 4 10^^ eV enregistres par AG AS A. Les carres 
rouges et les ronds verts representent les rayons cosmiques ayant des 
energies > ICP^ et (4-10) 10^^ eV respectivement. Les cercles plus grands 
indiquent les multiples (triplets en mauve). Ces evenements ont ete en- 
registres par le reseau d'AGASA de 20 km^ entre 17 fevrier 1990 et le 
31 juillet 2001, les angles zenithaux sont inferieurs a 45°. 
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1.5 Sources et mecanismes d'acceleration des RC 



La densite d'energie totale des RC mesures au-dessous de I'atmosphere est domi- 
nee par les particules avec des energies entre 1 et 10 GeV. A des energies au-dessous 
de ~ IGeV les intensites sont temporellement correlees avec I'activite solaire, ce 
qui est une evidence directe pour une origine solaire. A des hautes energies le flux 
observe montre une anti-correlation temporelle avec une activite solaire, indiquant 
une origine a I'exterieur du systeme solaire. Plusieurs arguments tels que la compo- 
sition, la production des rayons cosmiques secondaires suggerent que la quantite de 
RC entre IGeV et au moins jusqu'a la region du genou est confinees dans la galaxie 
et qui est probablement produite dans les restes de supernova (SNRs). Entre le 
genou et la cheville la situation est moins claire. Finalement, au dela de ~ lOEeV, 
les RC sont generalement supposes etre d'origine extragalactique, ceci est du a leur 
isotropic apparente. 

1.5.1 Mecanismes d'acceleration des RC charges 

De nombreux modeles tentent d'expliquer I'observation des rayons cosmiques 
au dela de 10^^ eV en introduisant des mecanismes d'acceleration qui permettent 
aux particules d'atteindre des energies tres elevees. 

La forme lisse du spectre evoquee plus haut laisse penser qu'il existe un nombre 
reduit de types de sources. II existe aujourd'hui un scenario assez communement 
accepte pour rendre compte des caracteristiques du rayonnement cosmique jus- 
qu'en dessous du genou. Les supernovae jouent ici un role preponderant. Celles-ci 
representent un moyen simple d'ejecter des particules relativistes dans le milieu 
interstellaire. Les abondances relatives dans le rayonnement cosmique correspon- 
draient ainsi de maniere assez naturelle aux abondances stellaires. 
Maintenant que Ton dispose d'une source, il faut accelerer ces particules de maniere 
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a obtenir le spectre observe. E. Fermi proposa en 1949 [H] un mecanisme ou les 
particules diffusees de maniere stochastique sur les irregularites magnetiques d'un 
nuage de gaz en mouvement a une vitesse V par rapport a la source gagnent sta- 
tistiquement une energie proportionnelle a (V"/c)^ a chaque collision ou c est 
la vitesse de la lumiere. C'est ce qu'on appelle le mecanisme de Fermi du second 
ordre. L'interet est que cela conduit naturellement a un spectre differentiel en loi de 
puissance. Etant donne que les vitesses tres faibles des nuages de gaz interstellaire 
et la faible densite de ceux-ci, I'efficacite du mecanisme est insuffisante. Quand on 
considere le meme phenomene, mais dans des ondes de choc, on obtient un gain 
d'energie en V/c on V est maintenant la vitesse du front de I'onde de choc. C'est 
le mecanisme de Fermi du premier ordre. Celui-ci conduit naturellement a une loi 
de puissance avec un indice spectral de -2. Ces ondes de choc apparaissent naturel- 
lement lorsque la matiere ejectee par une supernova rencontre le gaz environnant. 
De plus, les supernovae sont assez frequentes et energetiques pour maintenir la den- 
site locale d'energie contenue dans le rayonnement cosmique. Cette derniere pent 
etre obtenue a partir du flux observe et est approximativement pe = leV/cm^. Si 
Ton considere un temps de confinement des particules dans la galaxie r =~ 10^ ans 
on obtient une puissance d'injection necessaire : 



ou Vgaiaxie Gst Ic volumc dc la galaxie. Compte tenu de I'estimation de la puissance 
fournie par les ondes de choc de supernova, cela ne requiert qu'une efficacite de 
quelques pour-cents pour le mecanisme d'acceleration [I6j. 

Pour expliquer I'allure du spectre energetique a tres haute energie, le modele d'ac- 
celeration directe a ete propose. Dans ce cas, I'acceleration est rapide et due a 
I'existence d'un champ electromagnetique intense. Ainsi, un pulsar par exemple est 
le siege, a sa surface, de puissantes inductions magnetiques dues a la rotation rapide 
de la matiere condensee. Dans ce cas les particules peuvent etre produites et acce- 




(1.3) 
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champ 




Fig. 1.5 Emission d'un pulsar (modele du rotateur oblique). 

lerees grace a la force electromagnetique indmte(FEM). L'energie acquise depend 
du rayon de I'objet en rotation (R), du champ magnetique (B) a la surface et de la 
frequence de revolution (f). D'apres le modele du rotateur oblique [TTj, les pulsars 
sont interpretes comme etant des etoiles a neutrons en rotation autour d'un axe 
qui n'est pas aligne avec I'axe magnetique. Les particules chargees sont produites 
par extraction a la surface puis accelerees par la FEM qui se cree. A travers le 
rayonnement synchrotron elles emettent des 7. De cette maniere remission pent 
etre observee seulement lorsque le cone forme des deux axes est oriente vers la terre 
(figure (II. 5p . II faut signaler que ce mecanisme pent expliquer I'emission des rayons 
gamma ayant des energies de I'ordre du MeV ; par contre il ne pent expliquer celle 
des rayons gamma de haute energie. 

Ce mecanisme presente quelques problemes : d'une part, I'acceleration des parti- 
cules depend naturellement de la distance a laquelle I'influence du champ magne- 
tique cesse et d'autre part, il ne conduit pas d'une maniere claire a la forme du 
spectre energetique observe. 
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Accelerateurs cosmiques des particules 

De maniere generale pour etre accelerees aux plus hautes energies, les particules 
doivent rester confinees dans le site accelerateur pour pouvoir interagir avec celui-ci. 
Le confinement d'une particule dans un certain site d'acceleration va dependre de 
I'intensite du champ magnetique que celui-ci contient. En comparant le rayon de 
Larmor tl de cette particule a la taille L de I'objet astrophysique qui contient le 
champ B, il est possible d'estimer I'energie limite Ei jusqu'a laquelle la particule de 
charge Ze va pouvoir etre acceleree, avant de pouvoir s'echapper : 

Cette equation va nous permettre de dessiner un diagramme dit de Hillas [H], qui va 
representer les sites supposees capables d'accelerer des particules selon leur champ 
magnetique et de leur taille. La particule est acceleree tant qu'elle est confinee dans 
le site, la condition d'acceleration est done : 

B{nG)L{kpc) > E{EeV)/Z (1.5) 

Cela se traduit sur le diagramme (figure II. 6p par des sources potentiellement acce- 
leratrices qui doivent se trouver au-dessus des lignes representees. 
II faut done un compromis entre le champ magnetique qui doit etre suffisamment 
grand pour confiner les particules dans le site accelerateur, et la taille de ce dernier 
qui doit egalement etre suffisamment grand pour que les particules gagnent as- 
sez d'energie avant de s'echapper. Cela restreint deja les types d'objets candidats. 
Certains d'entre eux sont deja ecartes a cause de leur petite taille. Une exception 
a cette regie et un des candidats les plus serieux a I'acceleration des RC d'ultra 
haute energie sont les etoiles a neutrons (pulsars), qui compensent leur petite taille 
(quelques km) par des champs magnetiques tres intenses. Les restes de supernovae 
sont egalement exclus a cause de leur champ magnetique trop faible. Les objets 
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extragalactiques semblent les candidats les plus serieux, comme les noyaux actifs 
de galaxies, les lobes de radio galaxies et les sursauts gamma. Dans le cas ou les 
rayons de Larmor sont plus grands que la taille de la Galaxie, c'est une raison qui 
favorisent Thypothese de la nature extragalactique des rayons cosmiques au dela de 
la cheville. 

1.5.2 Sources des 7 de haute energie 

Les particules neutres telles que les 7 et de tres haute energie conservent la 
direction de la source et fournissent des informations precises sur la position de la 
source emettrice. Les neutrinos peuvent etre detectes par des experiences souter- 
raines grace a la detection des muons produits par les interactions de ce v avec la 
roche de couverture. Les 7 peuvent etre detectes par des detecteurs a la surface 
EAS "Extensive Air Shower" en etudiant les gerbes electromagnetiques en cascade. 
Les 7 de haute energie sont due probablement aux interactions des protons avec 
le milieu environnant d'une etoile a neutrons. Les protons acceleres, selon les mo- 
deles cites precedemment, jusqu'aux hautes energies peuvent interagir avec I'etoile 
accompagnatrice et le milieu interstellaire et donnent naissance aux ir" qui se desin- 
tegrent en deux 7. Des vr^ sont egalement produits et donnent apres desintegration 
des z/ et /i. 

Parmi les sources de 7 de haute energie nous citons les pulsars, qui emettent des 7 
dans la bande d'energie (l-lOO)TeV, dans la bande des rayons X nous trouvons le 
pulsar Crab associe a la nebuleuse Crab et Vela [20] . 

Les systemes binaires peuvent egalement etre responsables de remission des 7 de 
haute energie. Un systeme binaire est defini comme une etoile a neutrons en rota- 
tion autour de I'etoile massive dite "accompagnatrice". Parmi les sources qui ont ete 
signalees comme emettrices dans les bandes 7 et X nous citons Cyg-X3 et Hercule 
X-1. 
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Fig. 1.6 Diagramme de Hillas : les sites sont classes en fonction de leur 
taille et de leur champ magnetique. Les sources au-dessus des differentes 
lignes peuvent potentiellement accelerer les particules indiquees. 
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1.6 Interaction des RC avec I'atmosphere 



Des leur arrivee a la haute atmosphere, les RC interagissent avec les noyaux de 
ce dernier dormant naissance a un grand nombre de particules secondaires, essen- 
tiellement des mesons, qui peuvent re-interagir on se desintegrer dans I'atmosphere. 
En interagissant avec les noyaux de I'atmosphere, ces hadrons produisent d'autre 
particules en induisant une cascade formee de photons, d'electrons et de hadrons. 
Sur la figure (II. 7p . nous donnons une representation schematique des processus d'in- 
teraction d'un rayon cosmique primaire avec un noyau atomique dans I'atmosphere. 
Une cascade est formee : 

- d'une composante electromagnetique provenant de la desintegration rapide 
des 71° en photons ; chaque photon engendre a son tour une cascade electro- 
magnetique secondaire par I'intermediaire des eff'ets de creation de paires et 
de bremsstrahlung, pour chaque interaction du RC primaire charge, plus du 
tiers de I'energie est absorbe par la composante electromagnetique 

- d'une composante hadronique constituee principalement de pions et kaons 
charges qui se desintegrent pour atteindre le sol sous forme de muons et de 
neutrinos. 

Dans notre travail on s'interesse a la composante muonique. Cette derniere joue un 
role important dans I'etude des rayons cosmiques. Les muons out une petite section 
efficace d'interaction, ils interagissent faiblement et peuvent passer librement entre 
les noyaux atmospheriques et finissent par disparaitre en se desintegrant. Leur vie 
moyenne est de 2.2 x 10~^ s. En effet, les muons sont traditionnellement appeles la 
composante penetrante des rayons cosmiques, ils constituent la composante domi- 
nante des RC au niveau de la mer et ils peuvent etre detectes facilement par des 
detecteurs souterrains 

En effet, les protons et les neutrons interagissent fortement et ne peuvent traverser 
une longue distance dans I'atmosphere, les noyaux lourds fragmentent par collisions 
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Fig. 1.7 Representation schematique d'une gerbe produite par collision 
d'un rayon cosmique avec les noyaux atmospheriques. 
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avec les noyaux de Pair, les electrons et les photons subissent intensivement les pro- 
cessus electromagnetiques de perte d'energie et atteignent le niveau de la mer en 
nombre reduit ; enfin toutes les particules ayant une vie moyenne inferieure a 10~^*^ s 
se desintegrent dans la region de la haute atmosphere. Seuls les neutrinos arrivent en 
nombre eleve, mais a cause de leur faible section efficace d'interaction ils ne peuvent 
etre detectes aussi facilement que les muons par des detecteurs souterrains. 

1.6.1 Production des /i lors d'une cascade hadronique 

Les muons cosmiques sont produits a la suite de la desintegration des pions et 
des kaons charges qui out ete generes a leur tour par I'interaction des RC primaires 
avec les noyaux atmospheriques. La vie moyenne des pions est 2.6 x 10~* s, celle 
des kaons est 1.2 x 10~®s. Les principaux modes de desintegration et les rapports 
d' embranchement sont les suivants : 

+ i7^ 99.98770% 

jjT + 63% 
7r- + 7r° 2L% 

Les muons de basse energie sont generes durant la derniere phase de developpe- 
ment de la cascade. Ceux de haute energie, sont produits durant la phase initiale 
de la creation de la cascade, sont generalement tres penetrants et contiennent des 
informations sur les primaires qui les ont generes. 

Une etude cinematique simple dans le systeme du centre de masse (CM) puis 
dans le systeme du laboratoire (SL) de la desintegration tt — > permet d'obtenir 
I'expression de I'energie moyenne du muon dans le (SL) en fonction de I'energie du 
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vr, , et de son angle d'emission par rapport a la direction initiale du pion. 

= ^r/v/i^^ (1-6) 

avec Ef;^^ ^ 109.7 MeV. 
on encore 

= ^r^f = E^^^.E'^^/{m^c^) = 0.78 ■ (1.7) 

L'angle moyen par rapport a la direction du vr avec lequel est emis le muon dans le 
(SL) est de I'ordre de 10~^ radian (qui pent etre considere en premiere approxima- 
tion egal a 0). Ceci permet d'admettre que le muon conserve la direction du pion 
initial. Un raisonnement analogue pent etre fait egalement pour la desintegration 
des K. 

1.6.2 Production des /i lors d'une cascade electromagnetique 

Les 7 de haute energie interagissent avec I'atmosphere terrestre pour donner 
une cascade electromagnetique. Le processus de photoproduction pent egalement 
avoir lieu avec une faible probabilite selon : 

7 + noyau — >■ hadrons{Tr, K...) (1-8) 

Par la suite les vr et K se desintegrent pour donner des /i et des u. La section efficace 
de ce processus a ete mesuree par des experiences utilisant des accelerateurs jusqu'a 
des energies ~ 20 GeV des photons incidents. Au dela de la resonance, (0.1 - 0.8) 
GeV, cette section efficace est de I'ordre de 100 /x6/nucleon et augmente legerement 
pour des energies superieures a 10 GeV [23]. Cependant il a ete remarque que dans 
le cas de I'interaction 7p, I'extrapolation des donnees des sections efficaces aux 
faibles energies pent sous-estimer la probabilite de production des /i. 
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La probabilite qu'un 7 produise un hadron au lieu de la creation d'une paire 
e~^e~ a une energie ~ 20 GeV est donnee par le rapport des sections efficaces des 
deux processus [I3] : 

p^^ = ^j^h^ ~ 2.8 X 10-3 (1.9) 

La probabilite de production de hadron au lieu de la creation de paires est de I'ordre 
1/300. Les /i sont aussi crees par la production de paires /x^ mais en petite quantite. 
La section efficace de ce processus est donnee par : 

2 

m 

a^^^+f,- = Cr^^e+e-(— |) ~ 2 X 10 ^ ■ a^^e+e- (l-lO) 

Les cascades electromagnetiques riches en muons et qui out ete observees par cer- 
taines experiences peuvent trouver leur explication dans une grande section efficace 
de photoproduction, autrement il est difficile de donner une justification a ces re- 
sultats. 

II faut signaler enfin que pour avoir un /i de 1.4 TeV (energie minimale qu'un 
/i doit avoir pour atteindre MACRO) (voir ci-apres 1.6) le 7 primaire doit avoir 
environ 300 TeV, alors que pour une cascade hadronique il faut 20 TeV. 

1.7 Perte d'energie des muons dans la roche de 
MACRO. 

Les muons tres energetiques reussissent a traverser la roche de couverture sous 
laquelle est installe un detecteur comme MACRO. En traversant une epaisseur X, 
les /i perdent de I'energie d'une maniere continue ou discrete. 
Cette perte d'energie pent avoir lieu suivant 4 mecanismes possibles : ionisation, 
Bremsstrahlung, production directe de paires et interactions electromagnetique avec 
les noyaux. 



26 



La perte d'energie moyenne peut s'ecrire 



/dE\ _ UE\ ^ /dE\ ^ UE\ ^ UE\ 

La perte d'energie par ionisation domine aux faibles energies, alors que les autres 
mecanismes sont importants a des energies superieures au TeV. 

Dans le cas de I'ionisation, la perte d'energie est un processus continu qui peut 
etre decrit par I'equation : 

(^) " "^^'^ IniE^a./m^c')] (1.12) 

Emax est I'energie maximale transferee a un electron du milieu et qui est donnee 
par : 

El 

^'"'^^ " E, + O.Sm^cVme ^^'^^^ 
La relation (|1.12p peut etre approximee par 

on a ^ 2.0 MeV / g.cw? . 

En ce qui concerne la perte d'energie par Bremsstrahlung, la section efficace 
varie comme da/dv ~l/f ou v est la fraction d'energie cedee par le muon au photon 
durant le processus radiatif. La section efficace de creation de paires se maintient 
toujours finie et decroit rapidement, apres avoir atteint le maximum pour des v 
voisins de la limite cinematique Vmin = 4me/4ii^^, comme En negligeant les 

fluctuations pour ces deux processus discrets de la perte d'energie des /x, on obtient 
les courbes de perte d'energie correspondant aux divers processus et a I'ensemble 
(voir figure ( 11. 8p [T3]). On constate I'allure constante de la perte d'energie par 
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ionisation alors que les deux autres mecanismes croissent en fonction de I'energie 
selon une forme logarithmique. 

La perte d'energie totale est parametrisee par : 

dE . E 

{-^jtotale = (1-15) 

ax Xo 
avec Xo ~ 2.5 X 10^ g/cm"^ dans la roche. 

Une estimation de I'energie minimale qu'un fi doit avoir a la surface pour at- 
teindre une profondeur X est [T3] 



E^'"" = E^{e~o - 1) (1.16) 

Ec est appelee energie critique et correspond a la valeur pour laquelle les pertes 
d'energie par bremsstrahlung deviennent plus importantes que celles par ionisation. 
Les muons qui se propagent dans la roche ont Ec ~ SOOGeV^. Dans le cas de 
MACRO, pour traverser la roche de couverture (3500 m.w.e.J^, les muons doivent 
avoir une energie minimale de 1.4 TeV a la surface, ce qui correspond a des primaires 
d'energie 20 TeV. A cette profondeur le flux des /i atmospheriques se reduit d'un 
facteur 10^ par rapport a celui a la surface (c'est a dire 100 ix/cm^s). 



1.8 Technique de detection des rayons cosmiques 

Les methodes experimentales developpees par les physiciens depuis pres d'un 
siecle pour etudier le rayonnement cosmique sont conditionnees par le type de par- 
ticule que I'on cherche a detecter, mais surtout par la valeur du flux a I'energie 
consideree (voir flgure (11.91) ). Pour des energies inferieures a une centaine de TeV, 

^1 m.w.e (1 meter of water equivalent) est une notation utilisee pour exprimer la profondeur 
des laboratoires souterrains en metre d'eau equivalent ; (1 m.w.e= 1 hectogramme/ cm^) . 
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le flux est suffisamment eleve pour qu'un detection directe des particules soit pos- 
sible (~ 10 particules/m^/sr/s a 100 GeV). Au dessus de ratmosphere les detecteurs 
sont installes dans des ballons atmospherique, ou sur des satellites. Les detecteurs 
utilises sont soit passifs tels que les cliambres a emulsion, les films de rayons X et 
des detecteurs a trace nucleaires, soit actifs tels que : les calorimetres, les compteurs 
a scintillation et les compteurs Cherenkov. 

Pour les energies superieurs a quelques centaines de TeV le flux des RC est plus 
faible (inferieur a quelques particules/m^/sr/an). Ceci necessite une importante 
surface de detection, ce qui est problematique pour des experiences embarquees et 
exposees pour de breves periodes de temps (un a trois jours dans le cas des bal- 
lons). La detection directe des particules n'est plus possible et I'etude des RC ne 
pent se faire qu'a partir du sol en utilisant des detecteurs de grandes surfaces et 
qui sont exposes pour de longues periodes. Les experiences souterraines, comme 
MACRO, detectent la composante muonique (voir 12.111 . Sur les sommets des mon- 
tagnes sont places des experiences "EAS" qui etudient les cascades electromagne- 
tiques produites par des RC primaires d'energies variant de 10^"^ a 10^°eV. Parmi 
les detecteurs disposes au sol, nous citons I'experience Auger, c'est un detecteur 
Hybride utilisant des proprietes de la gerbes pour en analyser les caracteristiques 
et remonter au rayon cosmique initial. 
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Fig. 1.9 Schema indiquant le domaine d' application des differentes me- 
thodes de detection des rayons cosmiques. Un ordre de grandeur du flux 
est indique en fonction de Venergie. 
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Chapitre 2 



L 'experience MACRO 



2.1 Introduction 

Actuellement, la physique des accelerateurs nous permet d'atteindre des ener- 
gies de I'ordre de 14 TeV dans le centre de masse. Ce seuil est associe aux particules 
supersymetriques ou a des constituants des quarks et des leptons ; ceci pourra etre 
verifie dans I'avenir avec des accelerateurs puissants tel que le grand collisionneur 
de hadrons LHC (Large Hadron Collider) au CERN. Un deuxieme seuil pourrait 
etre associe a la grande unification qui se propose d'interpreter les trois interactions 
electromagnetique, nucleaire forte et faible et pouvant se manifester a des energies 
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de I'ordre de 10 eV. Parmi les consequences de ces theories, la desintegration du 
proton et I'existence de certaines particules exotiques tel que les monopoles ma- 
gnetiques, les nuclearites, etc.... L'observation de telles particules ne pent se faire 
aupres des accelerateurs. Dans ce cas, les experiences sans accelerateurs, en parti- 
culier les experiences souterraines, jouent un role tres important pour la recherche 
des particules traces et la violation de certaines lois de conservation. 
Les experiences souterraines sont souvent installees dans des sites miniers ou dans 
des tunnels, la roche couvrant I'experience joue le role de filtre qui absorbe les par- 
ticules chargees des RC dont le nombre est eleve au niveau du sol et peuvent rendre 
impossible la detection des evenements pen frequents. Les experiences souterraines 
sont nombreuses et sont destinees a differents sujets de recherche en physique fon- 
damentale. Neanmoins on pent les diviser en deux types, d'une part celles qui se 
basent sur I'etude des rayons cosmiques et d'autre part celles qui ne dependent pas 
de ces derniers. 

Les experiences souterraines qui ne sont pas destinees a I'etude des rayons cos- 
miques s'interessent generalement aux problemes de violations des nombres lepto- 
nique ou baryonique en particulier la desintegration du proton |24| . 

Les experiences souterraines basees sur les informations fournies par I'etude des 
rayons cosmiques out pour objectifs de resoudre les problemes relatifs : 

1. a I'astrophysique tel que : 

- I'etude de la composition et du spectre energetique des RC 

- I'etude de la variation temporelle et spatiale des RC, recherche des sources 
ponctuelles et detection des neutrinos de I'effondrement gravitationnel. 

- I'etude des neutrinos atmospheriques et leur oscillation. 

2. aux particules elementaires et a la recherche des particules exotiques telles 
que les monopoles magnetiques, les nuclearites, les Q-balles, les WIMPs, les 
particules supersymetriques, etc... 
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En general le nombre d'evenements attendus est tres petit ce qui necessite des 
detecteurs assez massifs et volumineux. 

Les detecteurs souterrains utilises sont : les scintillateurs liquides, les detecteurs 
Cherenkov a eau , les calorimetres avec reconstruction de traces, les tubes a streamer 
et les detecteurs plastiques. Pour I'etude des RC de haute energie (astronomic des 
muons et des neutrinos), on choisit generalement un des trois premiers types de 
detecteurs ou la combinaison de deux d'entre eux, afin d'augmenter la surface de 
detection et d'avoir une bonne resolution spatiale. Dans le cas de I'etude de la 
desintegration du proton, la competition se fait entre les detecteurs Cherenkov a 
eau et les calorimetres. 

Parmi les laboratoires souterrains nous citons : Le laboratoire Baksan en Rus- 
sie, I'experience Homestake, I'experience 1MB et Soudan 1-2 aux Etats Unis, I'ex- 
perience NUSEX sous le Mont Blac en France, I'experience Kolar Gold Field en 
Inde, I'experience Kamiokande et super-Kamiokande au japan enfin le laboratoire 
de Gran Sasso en italic. 

2.2 Generalites sur le Laboratoire National du Gran 
Sasso (LNGS) 

Le Laboratoire National du Gran Sasso est considere parmi les plus grands la- 
boratoires souterrains du monde, il est destine a la recherche en physique nucleaire, 
en physique des particules elementaires et en astrophysique. II comprend plus d'une 
dizaine d'experiences qui s'interessent a differents sujets de recherches [27] . 

Le LNGS est situe sous la chaine montagneuse du Gran Sasso au centre d'ltalie 
dans le tunnel autoroutier qui lie Rome a Teramo (voir figure (12. ip ). La latitude 
et la longitude du LNGS sont respectivement 42° 27' 09" Nord et 13° 34' 28" Est ; 



34 



<-L'Aquila 



Teramo-> 



Fig. 2.1 Vue generale du laboratoire souterrain. L'experience MACRO 
est localise dans le Hall B du lahoratoires Gran Sasso. 

I'altitude est d'environ 963m par rapport au niveau de la mer. 

La roche couvrant le laboratoire est constituee de calcaire avec une densite de 
2.71 ± 0.05(7/0771 , un numero atomique moyen 11.4 ± 0.2 et un poids atomique 
moyen de 22.9 ± 0.4 ; I'epaisseur moyenne de la roche est de I'ordre de 1400 m. Un 
muon arrivant au sommet de la montagne doit avoir au minimum une energie de 
1.4 TeV pour arriver au laboratoire souterrain [28] . 

Le LNGS contient trois grands Halls ayant chacun une dimension approximative 
de 100m x 17m x 17m [24]. 

2.3 Generalites sur le detecteur MACRO 

MACRO (Monopole Astrophysics and Cosmic Ray Observatory) [29] est situe 
au Laboratoire National du Gran Sasso (LNGS), il a ete considere comme etant 
un des plus grands detecteurs souterrains du monde [30]. II a ete realise en tenant 
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compte des objectifs exigeants : une acceptance suffisamment grande pour le flux 
de particules detectees, des techniques de detection varices (Scintillateur liquide, 
tube a streamer (TS) et le detecteur CR39) assurant une bonne resolution et une 
grande efficacite pour I'identification des diverses particules, et enfin un systeme 
d'acquisition performant. 

II a permis d'etudier difl'erents themes de recherches. 

MACRO est un detecteur a grande surface d'une longueur de 76.7 m, d'une 
largeur de 12 m et d'une hauteur de 9.3 m. II est constitue de six supermodules 
(SM) places I'un a cote de I'autre (figures ((221) et (ESD). Les six (SM) ont une 
structure identique. Un SM a une dimension de I'ordre de 12 m x 12 m x 9.3 m 
et il est forme de deux parties [3T] . 

Les deux modules sont separes par une distance de 30 cm qui est occupee en 
partie par le materiel supportant le detecteur. 

La partie inferieure est formee de deux plans horizontaux de scintillateurs li- 
quides separes entre eux de 4.8 m par dix plans de tubes a streamer limites. Entre 
les plans des tubes a streamer on trouve 60 cm d'absorbant, servant a eliminer le 
fond du au rayonnement 7 mous qui peuvent produire I'effet compton et de recon- 
naitre les muons ayant une energie > 1 GeV. Au milieu de chaque SM il y a un plan 
de detecteurs nucleaires solides a traces (CR39). Ce dernier est place egalement sur 
la partie laterale Est et frontale Nord. 

Afin d'augmenter I'acceptance du detecteur, les parois laterales et frontales sont 
constituees de six plans de tubes a streamer et un plan de scintillateurs. II faut noter 
que chaque supermodule pent fonctionner independamment des autres, permettant 
la continuite d'acquisition des donnees meme durant la reparation d'une partie du 
detecteur. 

La partie superieure de MACRO est appelee "attico" placee a une hauteur de 
4.8 m par rapport la partie inferieure. Elle est formee d'un plan de scintillateurs 
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horizontaux avec deux plans de tubes a streamer au dessus et deux autres au 
dessous. Comme pour la partie inferieure, L'attico a egalement des scintillateurs 
verticaux et des tubes a streamer lateraux. L'electronique et le systeme d'acquisition 
des donnees sont situes entre la partie inferieure et l'attico comme c'est indique sur 
la figure B . 

A I'interieur de MACRO est installe un detecteur a rayonnements de transition 
(TRD), ce module a pour objectif de mesurer I'energie des particules chargees qui 
le traversent. 



2.4 Les detecteurs constituant MACRO 

Le detecteur MACRO est constitue de divers subdetecteurs, (scintillateurs, 
tubes a streamer limites et le detecteur plastique a traces (CR39)). Dans ce pa- 
ragraphe, on va donner une description de ces detecteurs ainsi que leur mode de 
fonctionnement et leurs contributions a la detection des diff'erentes particules. 

2.4.1 Description des scintillateurs liquides 

Le role des scintillateurs dans I'experience MACRO consiste a mesurer la perte 
d'energie, le temps de vol, la direction et I'instant d'arrivee des particules en plus 
de la detection des neutrinos de I'effondrement gravitationnelle. 

MACRO contient 294 compteurs horizontaux et 182 compteurs verticaux. La 
masse totale du liquide scintillant est 600t. La partie inferieure du systeme de 
scintillateurs de MACRO de chaque SM contient 32 compteurs horizontaux et 21 
verticaux. Chaque compteur a une longueur de 12m dont 11m rempli d'un liquide 
scintillant transparent. 
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Fig. 2.2 Vue generale de I'experience MACRO. 



La figure (12. 4p donne une vue horizontale et verticale des scintillateurs horizontaux. 
Le scintillateur est constitue d'un recipient en PVC rempli de liquide scintillant, 
dont les parois internes sont enveloppees d'une couche de FEP-Tefion afin que ces 
parois deviennent totalement refiechissantes. Le liquide scintillant a la composition 
suivante [32] : 

- 96.4% d'huile minerale avec peu de parafine. 

- 3.6% de pseudocumene. 

- 1.44 mg/1 de bis-MSB. 

- 1.44 g/1 PPO 

La densite de ce melange est de 0.85g/cm^ et une longueur d'attenuation de ~12m. 
L'huile minerale est tres claire et a une longueur d'attenuation de 20 m. Le pseudo- 
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Fig. 2.3 Vite de Za section transversale du detecteur MACRO oil apparait 
tous les constituants du detecteur ainsi que leurs dispositions. 
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Fig. 2.4 Geometrie de Vune des extremites des compteurs a, scintillation 
horizontaux (gauche) et verticaux (droite). La partie centrale est remplie 
de liquide scintillant qui est vue par des PMT. 

cumene est un compose organique dont la composition chimique est (1,2,4-thrimethyl 
benzene) ; le PPO (2,5-diphenyl-oxazole) et le Bis-MSB (p-biso-methylstyrylbenzene) 

Les compteurs a scintillation horizontaux ont une dimension de 12 m x 25 cm 
X 75 cm chacun et ils sont mis sur trois plans horizontaux. La lumiere produite 
dans un compteur horizontal est recueillie aux extremites par un couple de photo- 
multiplicateurs (PMT) et les signaux issus de ces derniers sont additionnes. 

Les compteurs a scintillation verticaux ont une dimension de 12 mx25 cm x 50cm 
et couvrent les deux plans lateraux et le plan frontal de MACRO. Ces compteurs 
ne contiennent qu'un seul photomultiplicateur a chaque extremite. 

MACRO permet de reconstruire I'energie deposee dans un compteur avec une 
resolution as/E = 0.3/v^ (E est exprimee en MeV). La resolution temporelle 
est inferieure a Ins, ce qui permet de mesurer le temps de vol des particules et 
de determiner ensuite leur vitesse. On pent alors distinguer les muons provenant 
du haut de ceux provenant du bas et qui sont produits lors de I'interaction des 
neutrinos muonique avec la roche. 
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Les scintillateurs servent egalement a la detection des v^. produits durant I'ef- 
fondrement gravitationnel par I'intermediaire de la reaction 

z7g+p^e+ + ?i (2.1) 

dans le liquide scintillant. L'energie moyenne des antineutrinos de supernovae est 
(i?) ^ 12 MeV car la section efficace de la reaction 12.11 croit avec l'energie de 
neutrino. Cette reaction sera suivie, apres moderation du neutron dans le compteur, 
par la reaction de capture suivante : 

n + + 7 (2.2) 

avec = 2.2 MeV . Le temps de moderation des neutrons est de I'ordre de 10/is et 
le temps de capture est 180/is. Le positron donne naissance a un signal de particule 
rapide et permet a I'electronique de diminuer momentanement le seuil du meme 
compteur de 5MeV a IMeV pour un temps de 850/is afin de detecter le photon de 
2.2MeV^. Les deux signaux provenant du positron et du neutron sont les signatures 
caracteristique d'un tel evenement. 

2.4.2 Description des tubes a streamer limites 

Les compteurs tubes a streamer limites permettent de detecter les particules 
chargees rapides et lentes. Ces detecteurs sont aussi appeles detecteurs a recons- 
truction de traces. lis sont constitues de fils et de strips. Ces derniers sont inclines 
de 26.6° par rapport aux fils. Les fils permettent la mesure des coordonnees de 
la particule entrante dans le detecteur dans une direction appelee X, les strips a 
leur tour permettent la mesure des coordonnees de la particule dans une direction 
appelee D [33] (voir figure [275ll . 

Les tubes contiennent un melange de gaz constitue de 72% d'Helium et de 28% 
de n-pentane. Les proportions de cette composition permettent I'exploitation des 
eff'ets Drell et Penning pour la detection des MMs lents (/? = 10^^ — 10^'^). 
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Fig. 2.5 structure d'une chambre formee de 8 tubes a streamer avec 
les strips qui forment un angle de 26. par rapport aux fils. 



Les tubes a streamer limites sont des cellules ayant chacune une dimension de 
3 cm X 3 cm x 12 m, 8 cellules sont groupees dans une chambre en PVC ayant 
une epaisseur de 1.5 mm et une dimension de 25 cm x 3 cm x 12 m. L'anode 
est representee par un fil de 12 m de longueur et 100 /xm de diametre et qui sont 
lies aux extremites a un generateur de haute tension. La cathode est sous forme 
de graphite couvrant trois parois internes de la cellule et elle est liee a la terre. La 
tension d'alimentation typique des tubes a streamer limites est 4.2 kV. 

Pour chaque SM il y a 648 unites ou sont ranges 5184 fils. Les traces des par- 
ticules reconstruites sur les dix plans ont une resolution spatiale sur I'axe des X 
relatif aux fils de I'ordre de 1 cm, pour I'axe D la resolution relative aux strips est 
de I'ordre de 1.2 cm. On pent done tracer la trajectoire de la particule qui traverse 
les tubes a streamer avec une grande precision c'est pourquoi le systeme des tubes 
a streamer limite est appele "systeme de reconstruction de traces". 
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Fig. 2.6 disposition du detecteur track-etch. 

2.4.3 Description du detecteur nucleaire a traces de MA- 
CRO 

Le role principal du detecteur solide nucleaire a trace de MACRO est la recherche 
des monopoles magnetiques (MMs). II est place verticalement sur le front nord et 
sur la partie laterale Est et horizontalement au milieu de chaque SM inferieur. La 
surface totale occupee est de I'ordre de 1263 [33| . 

Ce detecteur se trouve range en unites de bases. Une unite de base est constituee 
d'un paquet de dimension 25 x 25 cm^ et formee de 3 feuilles de Lexan d'une epais- 
seur de 0.2 mm, 3 couches de CR39 d'epaisseur 1.4 mm et d'une feuille d'aluminium 
d'epaisseur 1 mm (voir figure (|2.6I) ). Les paquets sont mis sur des chariots mobiles 
en plastique afin de faciliter le retrait si un evenement est candidat a etre un MM 
par les scintillateurs et/ou les tubes a streamer limites. D'autre part le CR39 pent 
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etre utilise comme detecteur passif sans avoir recours aux triggers. 

C'est un detecteur plastique isolant, ayant pour formule chimique CH12H18O7, 
qui peut etre utilise dans plusieurs experiences en physique, en particulier la detec- 
tion des MMs recherche des particules ayant des charges fractionnelles, I'etude 
de la section efficace de fragmentation des ions relativistes [35] et la mesure de la 
concentration du radon |36| . 

La resolution en charge du detecteur CR39 est de I'ordre de a = 0.16e pour 
le balayage d'une seule surface. Pour dix mesures sur dix surfaces successives a = 
0.05e pour z = 6 et a = 0.07e pour z > 16. Cette resolution peut etre consideree 
comme etant la meilleure comparee a celles des autres detecteurs actifs. 

Le detecteur CR39 est sensible a la perte d'energie restreinte (REL) qui repre- 
sente la fraction d'energie concentree dans un cylindre de 0.01 /xm autour de la 
trace. La REL represente done une partie de la perte d'energie totale, I'autre partie 
d'energie est perdue sous forme de rayonnement 6 d'energie superieure a une va- 
leur specifique dependant du type du materiel utilise. Pour le detecteur CR39 cette 
energie specifique est de I'ordre de 200 eV. Le rayonnement S depose son energie 
loin de la trajectoire de la particule et ne peut done pas participer a la formation 
de la trace. 

Le mecanisme de formation de la trace latente depend de la particule nu- 
cleaire hautement ionisante. Un ion de numero atomique Z, en traversant un milieu 
condense se comporte comme une particule de charge effective positive [37] : 

Z,ff = Z^[l-exp(-130/3/Z2/3)] (2.3) 
si V est elevee, on a Z^jf ~ Z. 

Cela va produire une rupture au niveau de la chaine moleculaire et a la produc- 
tion de radicaux libres interagissant chimiquement. Les electrons de haute energie 
ne contribuent pas a la formation de la trace mais contribuent seulement a la perte 
d'energie dans un diametre de ~100 /im de la trajectoire de la particule. 
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2.4.4 Systeme d'acquisition des donnees de MACRO 

Le systeme d'acquisition des donnees de MACRO est un ensemble de calcula- 
teurs permettant la liaison a grande distance, la lecture des donnees relatives aux 
evenements ainsi que le controle du detecteur. 

L'acquisition est assuree par un reseau de six microvax et d'un vax central. Les 
microvax et le vax centrale sont lies par un reseau ETHERNET. Les microvax 
fonctionnent avec un systeme operatif VAXELN alors que le vax central qui est un 
VAX 4000/5000 fonctionne avec le systeme operatif VAX/VMS. 
Le systeme VAXELN optimise des operations temporelles et permet I'echange ra- 
pide et I'ecriture des messages entre les differents constituants du reseau. Les mi- 
crovax assurent le bon fonctionnement du detecteur, lisent et preselectionnent les 
evenements avant de les transmettre au vax central. lis sont branches en parallele 
avec les parties electroniques de I'appareillage a travers une interface de modules 
CAMAC et VME. lis s'occupent d'autre part des operations periodiques de calibra- 
tions des TS et sont affectes a certains controles du fonctionnement de I'appareil. 
Les trois premiers microvax controlent les informations des SM groupes en deux, 
les trois autres enregistrent les informations relatives aux neutrinos des supernova. 
Le systeme operatif VAX/ VMS qui gere le VAX central joue le role de serveur 
NETWORK et I/O pour le systeme operatif VAX/ VMS pour le systeme VAXELN 
qui gere les microvax et les composants du systeme VMS. II constitue ainsi le sup- 
port des histogrammes et rassemble les alarmes des differents microvax pour les 
memoriser sur des disques de haute capacite. II permet aussi la liaison du systeme 
d'acquisition avec I'exterieur par I'intermediaire de fibre optique. 



45 



2.4.5 Controle spontane d'acquisition 

Le systeme d'acquisition des donnees utilise par MACRO est represente sur la 
figure (12 .Tp . il est controle en temps reel par trois principaux processus s'activant 
sous le systeme VMS/VAX. 

- Un programme et un terminal jouant le role du console et a partir duquel 
seront executees les commandes de I'initialisation de I'acquisition. II permet 
le controle des diverses parties de I'appareil en temps reel (la frequence des 
triggers, les tensions d'alimentation, le temps mort...). On pent egalement 
voir les messages d'erreurs et les alarmes. 

- Un programme de presentation des histogrammes : "Histogram Presenter" 
permet de voir si toute les composantes de I'appareil fonctionnent normale- 
ment en visualisant des histogrammes de fonctionnement et compares a ceux 
de reference. Parmi les representations graphiques nous avons le temps de 
comptage de chaque trigger, le nombre de muons detectes par demi-lieure, le 
temps de vol et le temps mort de chaque microvax par demi-heure. 

- Un systeme de visualisation de la trajectoire des evenements dans les differents 
plans du detecteur "EVent Display" (EVD). Sur la figure (I2.8P on presente 
un groupe de muons visualise par EVD de MACRO. 

2.5 Les recherches en physique fondamentale de 
I'experience MACRO 

MACRO a ete destine a la recherche des monopoles magnetiques et des rayons 
cosmiques comme son nom I'indique. Cette recherche est tellement variee de telle 
sorte qu'il avait plus d'une dizaine de sujets qui ont ete traites par les physiciens de 
MACRO. Dans cette section on discute les differents sujets traites par le detecteur. 
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Fig. 2.7 Systeme d' acquisition des donnees utilise par MACRO. 
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Fig. 2.8 Groupe de 10 muons visualise par EVD du MACRO. 
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2.5.1 Detection des Monopoles magnetiques avec MACRO 

Dans ce paragraphe, on presente les derniers resultats de MACRO concernant 
I'acceptance et le flux des MMs |3^|39]. Pour cela il existe deux methodes d'ana- 
lyses : 

- I'analyse independante qui donne les resultats de recherches concernant chaque 
detecteur independamment de I'autre 

- I'analyse combinee concernant tout MACRO comme un seul detecteur. 

Ces etudes out ete menees pour difl'erentes vitesses des MMs. Sachant qu'apres une 
recherche de plusieurs annees aucun candidat MM n'a ete trouve, on presente les 
limites superieures sur le flux. 

A faibles vitesses, la recherche des MMs avec les scintillateurs ou les tubes a streamer 
est tres eflicace . Par contre elle est moins efiicace pour les MMs de grandes vitesses. 

Les resultats presentes ici sont relatives au MMs de charge unite. Les resultats 
sont valables en considerant une section efiicace de catalyse inferieure a 10 mb. On 
suppose que le flux des MMs est isotrope. 

a) Recherches des MMs avec les scintillateurs liquides : 

La recherche des MMs avec les scintillateurs liquides se fait sur une large gamme 
de vitesses ; ce qui necessite des methodes differentes pour la selection et I'analyse. 
On distingue trois domaines de vitesses pour la recherche des MMs : 

Pour les MMs de faibles vitesses, la perte d'energie des MMs pour 10~^ < j3 < 
5 X 10^^ est similaire a celle des protons ayant des vitesses analogues dans le liquide 
scintillant. La limite superieure sur le flux etablie pour des vitesses du MM pour 
10"^ < [3 < IQ-^ est de 3.4 x lO^^^ cm-'^ s"^ sr'^ [40]. 

Pour les vitesses moyennes des MMs, la lumiere produite par unite de longueur 
et le temps de vol pour parcourir MACRO sont utilises en meme temps. Dans 
cette procedure, le bruit de fond, principalement dii a la radioactivite naturelle et 
aux muons atmospheriques arretes dans le detecteur, est completement rejete. La 
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limite superieure sur le flux a 90% CL est 3.2 x 10 ^^cm ^sr ^ pour des vitesses 
1.2 X 10-3 <p < 10-\ 

Pour les MMs de grandes vitesses, la lumiere produite par perte d'energie des 
MMs de grande vitesse est collectee par les photomultiplicateurs. Pour ce type 
d'analyse le bruit de fond est dli essentiellement aux muons de haute energie et qui 
peuvent etre rejete puisque I'energie perdue par les muons dans les scintillateurs est 
nettement inferieure a celle que perdrait le MM. Cette analyse a permis de fixer une 
limite superieure sur le flux et qui est egale a 4 x 10~^^ cm^^s^^sr~^ pour /5 > 10^^. 

b) Recherches des MMs avec les tubes a streamer : 

La presence de I'helium dans les tubes a streamer donne une eflicacite de 100% vis 
a vis de la detection des MMs en exploitant I'efl'et Drell et Pening pour les faibles 
vitesses. L'analyse est basee sur la recherche de traces singulieres et la mesure de la 
vitesse du MM candidat . Seuls les plans horizontaux des tubes a streamer limites de 
la partie inferieure de MACRO ont ete utilises dans le trigger. Les plans de la partie 
superieure et les plans verticaux ont ete utilises pour la reconstruction des traces. La 
limite sur le flux obtenue par MACRO pour des vitesses 1.1 x 10~^ < /5 < 5 x 10"^ 
est 8.9 X 10-^^ cm-^s-^sr-^ [33j. 

c) Recherches des MMs avec le CR39 

Le CR39 est utilise dans MACRO comme un detecteur passif pour la detection des 
MMs, il pent etre egalement utilise pour la confirmation du passage d'un MM une 
fois sa detection est signalee par les scintillateurs et/ou les tubes a streamer limites. 

Le CR39 a ete soumis a des calibrations avec des ions lents et rapides afin de 
connaitre sa reponse lors du passage des particules. La limite superieure sur le 
flux trouvee est 6.8 x 10~^^cm~^s~^sr~^ pour des vitesses de I'ordre de /3 ~ 1 et 
10~^^cm~^s~^sr~^ pour des vitesses de I'ordre de /? ~ 10"'^. 

d) Analyse combinee 

L'analyse combinee est une analyse faite en utilisant les resultats des differents 
detecteurs constituant MACRO afin d'avoir une limite superieure combinee sur le 
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flux [n]. 

Cette analyse a ete faite pour 2774 jours avec une efiicacite moyenne de 77% et 
une acceptance de 3565 m'^sr ; la limite sur le flux trouvee est 7.6 x 10~^^ cm~'^s~^sr~^ 
pour des vitesses 5 x 10^^ < P < 0.99. Sur la figure (12. 9p on presente la limite supe- 
rieure sur le fiux des MMs du detecteur MACRO confrontee aux autres detecteurs. 
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Fig. 2.9 Limite superieure sur le flux des MMs a 90%CL avec g = go, 
obtenue par MACRO et par d'autre experiences. 
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2.5.2 Les Neutrinos des rayons cosmiques 

Les neutrinos sont caracterises par leur pouvoir de penetration et leur pro- 
pagation en ligne droite (sans subir les modifications dues au champ magnetique 
galactique). lis peuvent etre utilises comme moyen pour la recherche de sources as- 
trophysiques ponctuelles des rayons cosmiques. Ceci se fait en detectant les muons 
du bas issus de I'interaction des neutrinos avec la roche se trouvant au dessous de 
MACRO : 

i/^ + X ^ /i + y (2.4) 

MACRO pent egalement mesurer le flux des neutrinos du bas et le comparer avec 
les predictions des Monte Carlo afin de mettre en evidence I'existence d'un eventuel 
deficit de ces particules [41] . 

MACRO, grace aux scintillateurs, est egalement sensible aux antineutrinos des 
supernovae [42] [32] par I'intermediaire de leur interaction avec les protons libres se 
trouvant dans le liquide scintillant suivant la reaction suivante : 

X^+p — >n + e'^ (2.5) 

Le seuil energetique de detection atteint par MACRO est de I'ordre de 7 MeV. 

2.5.3 Les muons des rayons cosmiques 

La composante muonique est tres importante dans les recherches faites par MA- 
CRO. On distingue : 

- Temps d'arrivee des muons : I'analyse du temps d'arrivee des muons est 
une technique permettant d'etudier les differents mecanismes aux quels ont 
ete exposes les rayons cosmiques avant leur arrivee sur terre et les modulations 
introduites dans la distribution du temps d'arrivee de ces derniers [26] [43]. 
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- Astronomie des muons : cette analyse porte sur la recherche des sources as- 
trophysiques ponctuelles des rayons cosmiques, les resultats obtenus concernent 
la recherche de sources ponctuelles a travers tout le del (All Sky Survey) et 
I'etude de sources periodiques connues comme emettrices de rayons X telles 
que Cyg-X3 et Her -XI [26] [43]. 

- Variation saisonniere : il a ete montre que la composante secondaire des 
rayons cosmiques presente des variations en fonction du temps. Ces variations 
sont d'origine meteorologique [34] . 

- Ombre de la lune : cette analyse met en evidence un deficit des muons 
provenant de la direction de la lune. Elle permet de verifier la precision an- 
gulaire du detecteur MACRO. Les resultats obtenus en analysant 30 millions 
de muons montrent un deficit de I'ordre de 3.7a [44] . 

- Decoherence : c'est la distribution de la distance de separation de deux 
muons dans le cas des evenements multiples. Les donnees experimentales ont 
ete comparees avec le modele Monte Carlo HEMAS pour le flux des muons 
souterrains. Ce qui permet d'examiner les efl'ets systematiques possibles dus 
aux incertitudes sur la section efiicace d'interaction des primaires, la perte 
d'energie dans la roche et I'efl'et du champ geomagnetique [45] . 

- Decorrelation : la fonction de decorrelation des muons decrit la relation 
entre la separation angulaire et spatiale des muons souterrains. Comme pour 
la fonction de decoherence, elle est sensible aux previsions du modele d'in- 
teractions des rayons cosmiques de hautes energies et independante de la 
composition des primaires [46]. 




2.5.4 recherche des nuclearites et des WIMPS 

L'existence de nouvelle forme de matiere stable contenant un nombre egal de 
quarks up (u), down (d) et strange (s) a ete proposee par differents auteurs. Cette 
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matiere peut avoir des masses allant de quelques GeV a la masse d'une etoile a 
neutrons. De Rujula et Glashow ont suggere que cette matiere etrange peut exister 
dans les rayons cosmiques [79|. MACRO avec les detecteurs a scintillation et le 
CR39, peut etudier I'existence de cette matiere sur un large spectre de vitesses qui 
peut s'etendre de /3 = w/c = 1 jusqu'a la zone des /3 des nuclearites qui peuvent 
etre captures dans notre systeme solaire. 

La recherche des nuclearites a partir du detecteur MACRO peut se faire essentiel- 
lement par les scintillateurs et le detecteur CR39. Par contre, les tubes a streamer 
ne sont pas sensibles a ces particules du fait que le gaz se trouvant dans ces tubes 
a une faible densite. 

Les WIMPS (Weakly Interacting Massive Particles) sont des particules non 
baryoniques et representent un candidat pour la matiere obscure [48]. II font partie 
de la matiere galactique lumineuse qui a une contribution inferieure a 0.1% dans le 
densite de I'univers. 

Les WIMPS peuvent etre captures par le soleil ou par la terre et a travers les 
interactions avec les noyaux ils se concentrent au centre ou ils thermaliseront. Une 
paire de telles particules {w,w), peut s'annihiler en emettant des paires u v. Les 
neutrinos emis peuvent avoir une energie moyenne Ey ^ rriyu/S (quelques dizaines 
de GeV). La recherche de ces particules peut etre realisee en detectant les /i du bas 
provenant des neutrinos ayant la direction du soleil ou du centre de la terre. 
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Chapitre 3 



Analyse de la distribution des temps 
d'arrivee des muons de haute energie 



3.1 Introduction 

II est generalement admis que la distribution du temps d'arrivee des rayons 
cosmiques galactiques est un phenomene aleatoire. Cependant il existe quelques 
mecanismes introduisant des modulations dans cette distribution. La detection des 
rayons gamma de haute energie (10^'^ — 10^^ eV) provenant des sources galactiques 
ou extragalactiques a conduit aux developpement de differents modeles qui de- 
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crivent ces sources dans I'espace environnant (pulsars, systemes binaires, etc.). 
Weekes [49] a suggere que I'analyse des rayons cosmiques (RC) pent reveler la pre- 
sence de composantes dont I'origine est un pulsar ou un objet celeste similaire qui 
n'est pas entoure de nebulosite, ayant une courte periode et qui n'est pas situe a de 
grandes distances. De ce fait, une analyse des distributions temporelles des rayons 
cosmiques est importante. Dans le cas des rayons cosmiques primaires charges, les 
modulations des RC sont reduites ou eliminees par I'effet des champs magnetiques 
galactiques et extragalactiques presents entre la source et la terre. Ceci provoque 
la perte de la periodicite du signal provenant des sources situees a une distance 
inferieure a 150 pc. Les memes effets sont observes dans le cas des RC arrivant 
d'une supernovae. Les rayons cosmiques de haute energie gardent la periodicite du 
signal, ils peuvent preserver la modulation et done on ne pent exclure la presence 
de structures dans la distribution temporelle. II est done interessant d'examiner ces 
distributions independamment de toutes les hypotheses. 

Plusieurs recherches utilisant des detecteurs souterrains ou en surface exploi- 
tant "les gerbes electromagnetiques dans Pair (EAS)", n'ont pas reporte de modu- 
lation dans la distribution temporelle des RC de haute energie et par consequent 
confirment I'arrivee aleatoire des rayons cosmiques. Parmi ces etudes nous citons : 

- Morello et al. [50], a I'aide des mesures EAS avec deux differents seuils d'ener- 
gie E>50 GeV et E > 2 10^ GeV. 

- La collaboration NUSEX [51], est une experience souterraine sous le Mont 
blanc utilisant un detecteur compose de plans de Per separes par des plans de 
tube a streamer. Cette experience a mesure le flux des muons a une profondeur 
de 400 m.w.e. 

- La collaboration MACRO [26], a mesure le flux des muons singuliers et mul- 
tiples a I'aide du detecteur MACRO fonctionnant avec 2 et 4 supermodules 
seulement et en utilisant uniquement la partie inferieure du detecteur. 

^ 1 m.w .e{lmeterwaterequivalent) = Ihg/crn^ 
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Deux experiences ont confirme la presence d'une composante non aleatoire dans 
la distribution du temps d'arrivee des rayons cosmiques de haute energie : 

1. Bhat et al. [52] ont analyse les impulsions de la lumiere Cherenkov produites 
dans I'atmosphere par des RC d'energie primaire Ep > 100 TeV. lis ont ob- 
serve que la correlation se produit au temps t < AOs. L'experience a utilise 
deux photomultiplicateurs avec un grand angle d'observation a 2743 m d'al- 
titude. 

La figure (13 .11) donne la distribution des temps d'arrivee de deux evenements 
successifs et montre une deviation superieure a 5 a pour At < 40s. 

2. Badino et al. |53j, en utilisant un detecteur souterrain dans le tunnel du 
Mont blanc a une profondeur de 5000hg / cm^ de roche, ont observe un exces 
de muons en etudiant la distribution temporelle des muons. Les muons de- 
tectes ont une energie superieure k = 380Gey et sont originaires des RC 
primaires d'energie E^ ^ 50TeV. L'exces a ete observe pour un intervalle 
de temps separant deux evenements successifs d'environ 38s. Get exces de 
muons (1%) mis en evidence dans un intervalle de temps de quelques dizaines 
de secondes, revele la presence d'une composante modulee superposee a la 
distribution temporelle poissonnienne. 

La suite de ce chapitre, est consacre a I'etude des distributions de temps d'arrivee 
des muons cosmiques, realisee avec les donnees du detecteur MACRO fonctionnant 
avec I'ensemble de ces supermodules (six). Notre objectif est de conclure sur la 
nature de ces distributions et ainsi se situer par rapport aux resultats obtenus par 
les experiences citees ci-dessus. 

Un grand echantillon de muons a ete utilise compare a celui utilise dans les 
travaux ulterieurs de MACRO 1241 • La collecte de tel echantillon etait possible en 

r i ^ 

utilisant la periode durant laquelle le detecteur fonctionnait avec les six supermo- 
dules y compris la partie superieure du detecteur "attico". II faut noter que cette 
experience est en mesure de detecter environ 6 10^ muons/ans grace a sa grande 
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Fig. 3.1 Resultat de I'experience de Bhat et al. {5^ . L'histogramme 
montre la distribution observee des temps d'arrivee des evenements dans 
un bin de 10s. La courbe represente Vajustement par la fonction expo- 
nentiel theorique (en considerant que I'arrivee des evenements est com- 
pletement aleatoire). Une deviation de cette distribution a ete observee 
pour t < 40s. 
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acceptance. De ce fait, nous nous proposons d'etudier le temps d'arrivee des muons 
de differentes energies collectes durant la periode allant de 1995 jusqu'a juin 2000. 
Cette analyse porte sur trois types d'evenements : 

- evenements singuliers (multiplicite egale a 1). Un tel evenement correspond 

a un /i simple qui est reconstruit aussi bien sur les fils que sur les strips. Ces 

muons sont generes par des rayons cosmiques primaires ayant des energies de 

I'ordre de 20 TeV. Avec une telle energie, la particule chargee des RC n'est pas 

influencee par les modulations solaires. Ces muons nous permettent d'obtenir 

des informations sur les RC primaires ayant une energie dans le domaine 
1013 _ iqU gy 

- evenements doubles (multiplicite egale a 2). Un tel evenement correspond a 
deux /i paralleles et reconstruits sur les fils et sur les strips. De tels muons 
sont originaires surtout des primaires charges ayant une energie de I'ordre de 
200 TeV ou des 7 de haute energie mais avec une faible section eflficace. Cette 
zone d'energie est interessante puisque differents types de sources galactiques 
peuvent produire des rayons gamma de telles energies. 

- evenements multiples (multiplicite >3). C'est un groupe d'au moins 3 /i re- 
construits (sur les fils et les strips). lis correspondent aux primaires ayant une 
energie > 200 TeV. 

Dans notre etude nous proposons de comparer la distribution temporelle des 
muons a celle obtenue d'un processus poissonnien. 

3.2 Selection des donnees 

Vu I'importance de cette etude et sa grande sensibilite a toute anomalie dans 
la prise de donnees, ces dernieres ont ete selectionnees avec le maximum de soin en 
respectant des criteres de selection des runs. 

Le lot de donnees utilise correspond a la periode allant de Juin 1995 jusqu'a Mai 
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2000. Durant cette periode d'acquisition, les tubes a streamer ont presente un temps 
mort faible inferieure a 0.4%. 

Ce lot d'evenements est selectionne pour etudier les distributions temporelles des 
muons d'energie superieure a 1,4 TeV au sommet de la montagne, qui arrivent de 
la partie superieure du ciel et ceux qui proviennent des regions restreintes definies 
dans le systeme de laboratoire (cones). 



3.2.1 Criteres de selection 

Les criteres de selection auxquels les evenements et les runs doivent satisfaire 
sont : 

- La duree des runs doit etre superieure a 2 heures (la moyenne de la duree des 
runs est d'environ 5 heures). Ceux qui ont une duree inferieure sont en general 
des runs de calibration on des runs qui presentent des problemes (electronique, 
manque de courant, erreurs du systeme d'acquisition, etc.). 

- La frequence des muons R est choisie dans I'intervalle 840 < R < 960 n/h 
et situee dans un intervalle ±2cr autour de la valeur moyenne. Les runs qui 
presentent une frequence elevee ont ete analyses separement pour verifier s'il 
s'agit d'une fluctuation due au bruit electronique ou d'un eventuel "burst" de 
muons. Ce critere nous a permis d'eliminer les runs ou I'acquisition n'a pas 
ete realisee avec tons les 6 SM ou ceux qui presentent des problemes. 

- Un temps mort d'acquisition inferieure a 0.4% pour chaque microvax. Ce qui 
permet d'eliminer les runs qui presentent un exces de bruits provenant de 
I'electronique ou des decharges continues des tubes. 

- Le temps d'arrivee des muons a ete mesure avec une horloge atomique d'une 
precision d'environ absolue et O.S/xs relative et ne doit presenter aucune 
anomalie. 
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- Absence de probleme dans le systeme du gaz dans les tubes a streamer. Un 
mauvais fonctionnement du systeme du gaz pent avoir un effet sur I'efficacite 
des tubes et par consequent sur la frequence d'arrivee des /i detectes. 

- L'efficacite de detection des tubes a streamer doit etre superieure a 90% et 
86% pour les fils et les strips respectivement. 

La figure fl3.2[ a) montre le nombre de muons en fonction du nombre de hits utilises 
pour la reconstruction de la trace des muons singuliers. lis faut signaler que seuls 
les 10 premiers plans ont ete utilises. On constate que la majorite des traces out ete 
reconstruites avec plus de 8 hits, ce qui donne une idee sur la grande efficacite des 
tubes a streamer. Les muons reconstruit avec moins de 7 hits, sont generalement les 
muons qui traversent deux modules juxtaposes ou qui sortent des parois lateraux 
du detecteur. 

Sur la figure (I3.2[ b) est presente le nombre de run en fonction du nombre moyen 
de plans des tubes a streamer qui ont contribue a la reconstruction de la trace. On 
constate, la aussi, que le nombre moyen de plan est 9.4 qui est une valeur impor- 
tante confirmant la grande efficacite de detection des tubes. 

La figure (13.31) . donne l'efficacite des fils de chaque run selectionne (la valeur 
moyenne de cette efficacite est de I'ordre de 94%). Cette efficacite est obtenue en se- 
lectionnant les muons qui traversent uniquement les dix plans des tubes a streamer, 
sans tenir compte des muons qui sortent lateralement, et en divisant le nombre de 
points qui contribuent a la reconstruction des traces par le nombre d'evenements 
par run et par le nombre de plans. 

D'autres criteres de selection ont ete appliques, un evenement de muon singulier 
doit avoir une trace avec une multiplicite egale a 1 sur les fils et sur les strips ; un 
evenement de muon double doit avoir deux traces sur le fil et sur le strip et les 
muons multiples ont des multiplicites entre 3 et 6 et reconstruits sur les fils et sur 
les strips. 

Nous presentons dans la figure (13.41) la frequence des muons singuliers en fonction 
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Fig. 3.2 (a) Nombre de muons en fonction du nombre de points utilises 
pour la reconstruction de la trace ; (b ) Nombre moyen de plans ayant 
contribues a la formation de la trace. Seuls les 10 premiers plans ont 
ete utilises pour cette etude. 
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Fig. 3.3 L'efficacite des fils des tubes a streamer pour chaque run selec- 
tionne. Elle est de I'ordre de 94%- 
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du numero des runs avant qu'aucune selection n'a ete appliquee (a) et apres les 
selections (b). 

Sur la figure f l3.4[ b) on constate la presence d'une variation reguliere et qui est la 
signature d'une variation saisonniere dans le flux de muons due aux changements 
de la temperature et de la densite de ratmosphere. En efl'et les muons etudies sont 
produits par la desintegration des mesons n et k dans I'atmosphere. Ainsi, le flux de 
muons depend du rapport entre la probabilite de desintegration et celle de I'interac- 
tion des mesons parents avec I'atmosphere et qui sont sensibles a la densite et a la 
temperature de I'atmosphere. Le flux de muons diminue en hiver, ou la temperature 
est minimale et I'atmosphere est plus dense, et augmente en ete. L'amplitude de 
cette variation est de I'ordre de ±2% et elle est presentee sur la figure (13 .Sp [54] . 

L'etude des variations journaliere du temps solaire et sideral du flux des particules 
montre que le detecteur presente une haute sensibilite pour l'etude des tres petites 
variations du flux de muons (de I'ordre de 10~^) ^55]. Elle montre bien la grande 
sensibilite de MACRO meme aux faibles variations du flux de muons, ce qui est 
important pour l'etude des variations des temps d'arrivee des muons et pour la 
recherche des sources astrophysiques des RC. 

Dans la figure (13. Gi b) nous avons presente le nombre de runs en fonction de la 
frequence d'arrivee de muons pour chaque run. La frequence moyenne est egale a 
897.7 

La duree des runs est donnee sur la figure (13. 61 a). La duree moyenne est de I'ordre 
de 5 heures. 
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Fig. 3.4 Nombre de muons par heure en fonction du numero des runs ; 
(a) avant la selection; (b) apres la selection. On constate bien la presence 
de la variation saisonniere dans le flux de muons. 
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Fig. 3.5 Variation saisonniere du flux de muons [25]. On observe une 
diminution du flux de muons en hiver et une augmentation en ete. L 'am- 
plitude de cette variation est de I'ordre de ±2% 
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Fig. 3.6 (a) Duree des runs selectionnes. (b) Nombre de runs en fonc- 
tion de la frequence d'arrivee des /i. La frequence moyenne est 898ii/h. 
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Fig. 3.7 Distribution des temps d'arrivee entre deux muons singuliers 
consecutifs pour dti < 1 s, notons I'effet du temps mort sur les premiers 
bins. 

3.3 EfFet du temps mort 

Apres chaque trigger il existe un temps mort de I'ordre de 100 ms lorsque les 
scintillateurs etaient en acquisition. L'effet de ce temps mort est apparent dans le 
cas des distributions (ti — to) pour (ti — to) < 100 ms, figure (13. 7p . L'etude de I'eff'et 
du temps mort realisee sur un echantillon de 4 10^ muons detectes par MACRO 
|26| . montre que I'ajustement avec ou sans correction donne le meme resultat. 
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3.4 Analyse des correlations temporelle des muons 



3.4.1 Generalites 

Comme il a ete deja mentionne, il est prevu que les rayons cosmiques ont un 
temps d'arrivee aleatoire a cause des deflexions des particules chargees des rayons 
cosmiques par les champs magnetiques interstellaires, cependant il se pent que cer- 
tains mecanismes introduisent des correlations dans le temps d'arrivee, par exemple 
un rayon cosmique originaire d'un pulsar proche pent introduire une modulation 
aux rayons cosmiques avec une frequence egale a la frequence de rotation du pul- 
sar [49]. Des effets similaires peuvent apparaitre avec d'autres types de sources de 
rayons cosmiques. 

Pour chaque muon qui arrive au temps t^, nous analysons la difference de son 
temps d'arrivee avec le second muon et les quatre muons consecutifs : ti — toi^dti), 
t2—to{dt2), t^—to^dts), t4—to{dt4^), t^—toi^dt^) comme c'est deflnit sur la flgure (I3.8p . 
Cette analyse est effectuee en comparant la distribution temporelle de ces muons 
a une distribution poissonnienne. On pent prevoir I'observation des clusters dans 
le temps des evenements generes par des particules chargees ou neutres produites 
dans une emission periodique d'une source proche. 

Le flux de rayons cosmiques, en I'absence de modulations dans les mecanismes de 
production et de propagation, est distribue d'une maniere aleatoire dans le temps. 
Par consequent nous comparons la distribution des retards entre un fj, et le n**^"^"^ /x 
(n=l,2,...) successif avec la fonction Gamma qui caracterise les phenomenes alea- 
toires. 

La fonction Gamma d'ordre M est deflnie comme suit [5^ : 

G{x;N,X,M) = NX ^ (M-1)! ' ^<^<^ (3-1) 
Ou A et M sont deux constantes reelles positives et N un facteur de normalisation. 
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Fig. 3.8 Definition des quantites dti 
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Pour des valeurs entieres de M, la fonction Gamma est appelee Erlangienne d'ordre 
M et elle est equivalente a un ensemble de distributions exponentielles. 

Pour M=l, I'equation (13.ip se reduit a la fonction exponentielle simple : 

A^, A, 1) = iVAe-^^ (3.2) 

et sera appliquee a la distribution de la difference entre les temps d'arrivee d'un ii 
et du suivant immediat {dti). 

Une caracteristique de la fonction Gamma est que, au contraire de M, le parametre 
A doit rester constant. Dans nos calculs, nous avons determine A dans le cas ou 
M=l. Pour les M>1, nous avons laisse le parametre A libre de varier. 
Un meilleur ajustement avec la fonction Gamma permet d'obtenir le coefficient 
/DoF voisin de 1. L'ajustement est effectue en utilisant le programme de mini- 
misation MINUIT de la librairie CERN [57|. 

3.4.2 Analyse des Muons arrivant de toutes les directions 

Nous avons effectue I'analyse sur un large lot de donnees constitue de 8.6 10^ 
de muons simples, 0.46 10^ muons doubles et 0.08 10^ muons multiples qui arrivent 
de la partie superieure du ciel. Cette partie est definie en coordonnees celestes 
par la bande limitee en ascension droite par 0° < a < 360° et en declinaison par 
-20° < (5 < 90° [24]. 

Nous presentons sur la figure (13.91) la distribution temporelle separants deux muons 
singuliers consecutifs (dti). L'ajustement des donnees experimentales donne une 
valeur de y^jDoF ^ 1, les resultats sont repertories dans le tableau (13.11) . Ceci 
indique une adaptation adequate avec la fonction theorique prevue. Par consequent 
la distribution des intervalles des temps d'arrivee entre deux muons consecutifs est 
compatible avec I'hypothese qui exclut la presence de toute correlation temporelle. 
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Fig. 3.9 Distribution du temps d'arrivee des /i singuliers (points). La 
courbe en pointilles est le resultat de I'ajustement des donnees par la 
fonction Gamma d'ordre 1. On constate le bon accord entre les don- 
nees experimentales et I'ajustement avec la fonction Gamma d'ordre 1, 
ce qui montre une arrivee aleatoire dans le temps d'arrivee des muons 
cosmiques. 
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Fig. 3.10 Distribution du temps d'arrivee des ji doubles. La courbe en 
pointilles est le resultat de I'ajustement des donnees par la fonction 
Gamma d'ordre 1. 
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Fig. 3.11 Distribution du temps d'arrivee des /i multiples. La courbe 
en pointilles est le resultat de I'ajustement des donnees par la fonction 
Gamma d'ordre 1. 
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L'absence d'une eventuelle correlation temporelle sur dti, nous a pousse a etu- 
dier les ordres superieurs tels que dt2, dt^, dt^ et dt^. 

Les resultats relatifs aux muons singuliers, doubles et multiples sont presentes sur 
les figures (13.121) . f l3.13p et (13.14p respectivement. Pour chaque ordre, I'ajustement 
a ete eff'ectue avec la fonction gamma d'ordre correspondant (M=2 pour dt2, M=3 
pour dts, etc.). Notons que I'ajustement s'accorde parfaitement avec les donnees 
experimentales, ce qui nous permet de deduire qu'aucune correlation temporelle n'a 
ete observee et de confirmer les resultats obtenues ulterieurement par la collabora- 
tion MACRO [26]. 

3.4.3 Analyse des muons arrivant des cones 

Les detecteurs de grande surface fixes par rapport a la terre nous permettent 
d'observer des petites anisotropics. Avec sa grande acceptance, MACRO nous per- 
met d'observer des directions qui changent avec la rotation de la terre et done une 
bande limitee et bien definie dans I'espace sera decrite durant les 24 heures side- 
rales. 

Pour chercher I'existence d'une composante non aleatoire, trois zones ayant des 
declinaisons separees dans le temps out ete choisies. Ainsi nous avons mesure les 
distributions des temps qui separent I'arrivee des /x dans trois fenetres spatiales 
definies en coordonnees locales et centrees sur des valeurs diff'erentes de Tangle azi- 
mutale. De ce fait, les trois bandes auront des declinaisons differentes, non correlees 
et se suivent dans le temps. 

Nous repetons la meme analyse pour I'etude des correlations des muons qui ar- 
rivent par des directions definies par les angles zenithal et azimutale. Dans ce cas 
le detecteur joue le role d'un telescope pour observer differentes region du ciel. Les 
cones (1) et (2) choisis coincident avec le maximum d'intensite de la distribution 
angulaire des muons en zenith et azimut ; le cone (3) a ete choisi pour couvrir la 
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Fig. 3.12 Distribution du temps d'arrivee des ji singuliers. La courbe 
en pointilles est le resultat de I'ajustement des donnees par la fonction 
Gamma d'ordre 2,3,4,5. 
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Fig. 3.13 Distribution du temps d'arrivee des ji doubles. La courbe en 
pointilles est le resultat de I'ajustement des donnees par la fonction 
Gamma d'ordre 2,3,4,5. 
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Fig. 3.14 Distribution du temps d'arrivee des ji multiples. La courbe 
en pointilles est le resultat de I'ajustement des donnees par la fonction 
Gamma d'ordre 2,3,4,5. 
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region de declinaison centree sur la direction de Cygnus-X3. Dans la figure (13.151) 
nous presentons les distributions des muons en fonction des angles azimutal (a) et 
zenithal (b) pour les muons singuliers. 
Pour le premier cone clioisi avec 25 < zenith < 45 et 20 < azimut < 40, le 
nombre total des muons singuliers selectionnes est 280 10^. Pour le deuxieme cone, 
25 < zenith < 45 et 140 < azimut < 160, nous avons selectionne 350 10^ fi. Pour 
le troisieme cone, nous avons choisi 60 < azimut < 100 et 260 < azimut < 300 
et 30 < zenith < 50 et le nombre total des muons singuliers dans ce cas est de 
520 10^. Les declinaisons et les ascensions droites des muons singuliers arrivant des 
cones out ete calculees en utilisant les coordonnees locales et le temps d'arrivee des 
evenements. Dans la figure (13 .16^ nous presentons les bandes de declinaison des 
muons singuliers selectionnes pour les diff'erents cones, les cercles indiques le plan 
galactique. 

Dans la figure ( 13.17^ nous montrons les distributions des temps d'arrivee des muons 
singuliers venant du cone 1 (cercles pleins), cone 2 (triangles) et du cone 3 (cercles 
vides). 

Les parametres d'ajustements des distributions experimentales ti — to avec la fonc- 
tion Gamma d'ordre 1 (Eq. I3.2p sont presentes sur le tableau (13.11) . On constate 
un bon accord entre les distributions theoriques et experimentales, ce qui exclut la 
presence de toute sorte de modulation. 
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Fig. 3.15 Distributions des angles Zenith (h) et Azimut (a) des muons 
singuliers. 
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Fig. 3.16 Distribution des muons singuliers venant des cones selection- 
nes en fonction de la declinaison et de I'ascension droite. Les cercles 
ouverts indiquent le plan galactique. 
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Fig. 3.17 Distribution du temps d'arrivee des jji singuliers consecutifs 
venant du cone 1 (cercles pleins), cone 2 (triangles) et du cone 3 (cercles 
vides). Les courbes representent le resultat de I'ajustement des donnees 
par la fonction Gamma d'ordre 1. 
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1.36 
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180.2±0.3 


196.0±0.4 


0.93±0.07 


0.98 


Cone 1 


276.3 ±0.6 


121.58±0.03 


1.003±0.004 


1.02 


Cone 2 


353±1 


95.9±0.1 


1.003±0.006 


0.82 


Cone 3 


522±1 


63.82±0.3 


1.004 ±0.004 


0.91 



Tab. 3.1 Parametres d'ajustements des distributions experimentales 
(ti — to) avec la fonction Gamma d'ordre 1 pour les fi singuliers, doubles 
et multiples arrivant de tout le del et des 3 cones. 

3.5 Test de Kolomogorov-Smirnov 

Afin de chercher des structures possibles dans les distributions des temps d'ar- 
rivee, nous avons aussi utilise le test de Kolomogorov-Smirnov [59]. Ce test nous 
permet de tester Thypothese Hq selon laquelle les donnees observees sont engendrees 
par une loi de probabilite theorique consideree comme etant un modele convenable. 
Le test compare la distribution cumulative F{x) des donnees experimentales avec 
la distribution theorique aleatoire H{x). 

La mesure de la deviation est d = max \ H{x) — F{x)\, ou F{x) et H{x) sont les 
distributions cumulatives de f{x) (donnees) et h{x) (attendue) respectivement. En 
terme de qualite du test, F{x) s'accorde avec H{x,X), ou A est prise des donnees, 
avec une probabilite de compatibilite, entre les distributions mesurees et attendues, 
donnee par : 

Pkid > observee) = QksiVlN)d) (3.3) 
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ou : 

+00 

Qfc,(x) = 2 5^(-irie-^V (3.4) 
i=i 

La probabilite des tests pour les distributions {ti — to) est donnee dans le tableau 
(13 ■21) . Ces resultats sont en accord avec les distributions aleatoires, bien que dans 
le cas du cone 3, un disaccord est remarque (a un niveau de 1 a) et pent produire 
une faible probabilite. Ceci pent etre du a une possible augmentation du rapport 
d'evenements pour (ti — to) = 500 s (voir figure (13.171) ). mais les faibles statistiques 
ne peuvent conduire a une conclusion assez claire. 



Selection 


Pr. K-S 


singuliers 


0.99 


/i doubles 


0.95 


H multiples 


0.99 


Cone 1 


0.99 


Cone 2 


0.77 


Cone 3 


0.38 



Tab. 3.2 Probabilite de Kolmogorov-Smirnov pour les fi singuliers, 
doubles et multiples arrivant de tout le ciel et des 3 cones. 
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3.6 Conclusion 



Nous avons presente les resultats de I'etude des distributions des muons cos- 
miques avec une energie plus grande que 1.3 TeV au sommet de la montagne de 
Gran Sasso. Les donnees analysees ont ete collectees a I'aide du systeme des tubes 
a streamer du detecteur MACRO dans sa configuration complete. 
Les muons singuliers, doubles et multiples arrivant de toutes les directions ainsi que 
ceux venant des zones selectionnees formant des cones ont ete consideres [60j. 
Les resultats de notre etude concernant les differentes distributions temporelles 
montrent que les donnees experimentales sont compatibles avec la fonction Gamma 
d'ordre M. Ce qui nous a permis d'affirmer que nos resultats sont en accord avec 
I'arrivee aleatoire des muons dans notre detecteur et permet ainsi d'exclure la pre- 
sence de toute composante, qui ne soit pas aleatoire dans les temps d'arrivee des fi, 
similaire a celle observee par Bath et al. [52], ou par Badino et al. [53], et ceci pour 
des energies de I'ordre de 20 TeV (// singuliers) et superieures a 20TeV (/i doubles et 
multiples). L'etude des distributions temporelles des fx arrivant des cones, a abouti 
aux memes conclusions. 

Le test de Kolmogorov-Smirnov, nous a permis de mesurer la deviation entre la 
distribution theorique et celle experimentale. Aucune deviation n'a ete mise en evi- 
dence dans notre echantillon de donnees, sauf dans le cas du cone 3 et ceci est dii 
a la faible statistique. 

Dans cette etude nous n'avons mis en evidence I'existence d'aucun signal qui pent 
etre la signature d'une source discrete tels que les pulsars. 
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Chapitre 4 



Recherche des variations du flux de 



muons 



4.0.1 Introduction 

Durant presque dix ans de prise de donnees a une profondeur de 3800 m.w.e, 
MACRO a rassemble un des grands echantillons de muons coUectes par les expe- 
riences souterraines de ce type. Le flux de muons pent presenter des variations 
d'origine galactique, solaire ou meme terrestre. Lorsque ces variations ne sont pas 
aleatoires, elles peuvent reveler la presence des modulations an niveau du signal 
devoilant ainsi I'existence des objets galactiques emetteurs. Les clusters de muons 
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souterrains peuvent etre produites par des evenements violents comme les "Gamma 
Ray Burst" sur I'echelle de temps de quelques secondes, ou par des variations me- 
teorologiques soudaines sur I'echelle de temps de plusieurs heures. Les modulations 
periodiques sont lies a la variation de la temperature de la haute atmosphere, qui 
diminue durant la nuit (hiver) produisant les variations journalieres (saisonnieres) 
de la densite de I'atmosphere et cependant du flux de muons. 
Le temps d'arrivee des muons cosmiques de haute energie a ete analyse precedem- 
ment, ou nous avons pu confirmer I'absence de modulations dans le flux de muons. 
Cette etude va etre argumentee dans ce chapitre par la recherche de possibles va- 
riations de flux. Ainsi deux methodes sont utilisees : 

- la recherche d'eventuel cluster d'evenements (groupement d'evenements) dans 
le flux de muons. 

- la recherche de variations periodiques. 



4.1 Recherche de cluster d'evenements 

On se propose de chercher la presence des clusters d'evenements dans le flux de 
muons. L'analyse de ces clusters nous permet d'estimer sa signiflcation statistique, 
ce qui nous permet de conclure sur sa nature qui pent etre un resultat du hasard ou 
a I'observation d'un exces d'evenements originaire d'un objet emetteur des rayons 
cosmiques. Dans ce dernier cas, qui presente generalement beaucoup d'interet, il 
est important de trouver et d'expliquer les facteurs qui out conduit a sa realisation. 
La presence de groupement d'evenements dans le flux de muons cosmiques, indique 
I'existence d'une composante emettrice des rayons cosmiques dans une direction 
privilegiee. 

Les perturbations locales dans une distribution theorique d'une variable donnee, 
peuvent signaler la presence d'une modulation, ou d'une composante en disaccord 
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avec le modele fondamental, utilise pour decrire les donnees experimentales (I'hypo- 
these nulle). En physique de hautes energies, ce phenomene est relie generalement a 
I'apparition d'une resonance inattendue, et/ou a la determination de la signification 
statistique de I'exces. Les anomalies globales sont frequemment traitees a I'aide du 
test de Kolmogorov-Smirnov et ses extensions. Cependant, le pouvoir de ce test 
est considerablement reduit dans le cas des perturbations locales. Inversement, le 
test utilise un binage de I'intervalle et compare le contenu de chaque bin avec 
celui theorique sous I'hypothese nulle. Ce test est mieux adapte aux perturbations 
locales. Notons que le test est adapte seulement dans le cas d'un binnage fixe a 
priori. En general, les perturbations locales partagees entre differents bins, sont 
moins marquees que ceux oil le cluster d'evenements est situe dans un seul bin. 
De ce fait, les techniques de scan avec une fenetre de longueur fixe semble les plus 
appropriees. 

Ces techniques sont generalement utilisees dans plusieurs domaines de recherche 
tels que la bio-informatique, la medecine [62], la physique des particules [69] ainsi 
que dans des applications en astrophysique [63J. 



4.1.1 Methode scan statistics 

La methode de balayage scan statistics est I'une des plus puissantes methodes 
utilisees pour analyser I'apparition de cluster d'evenements. Elle presente un moyen 
tres utile pour signaler la presence de perturbation dans le modele de probabilite 
fondamentale qui decrit les donnees experimentales. Elle est basee sur I'echantillon- 
nage libre [61] . 

Soit X une variable continue de I'intervalle [A, B] et qui obeit a un processus de pois- 
son, notons A la valeur moyenne par unite d'intervalle. La probabilite de trouver 
Yx{w) evenements dans un intervalle [x,x + w] est 

Prob{Y,{Lo) = k) = e-"--^ A; = 0, 1, 2, ... (4.1) 
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Le nombre d'evenements dans les intervalles disjoints est independamment distri- 
bue. Nous appelons scan statistics (SS) le plus grand nombre d'evenements trouve 
dans un sous-intervalle de [A, B] de longueur u : 

S{u;)^ max {Y^{w)} (4.2) 

A<x<iB—uj 

La probabilite pour que le nombre d'evenements dans une fenetre balayee n'atteint 
jamais k sera donnee, selon [62], par : 



Ou A = B — A et\e suffixe "*" indique que les probabilites non conditionnelles sont 
considerees, i.e. que le nombre global d'evenements N dans I'intervalle fluctue selon 
I'Eq. f l4.ll) avec = A. La forme exacte de I'Eq. (14. 3p pent etre exprimee sous forme 
d'une somme des produits de deux determinants [64] [M]- La sommation se fait sur 
un ensemble V constitue par 2H + 1 partitions de N dont les elements mj de chaque 
partition sont des entiers non-negatifs satisfaisants la condition mj + mj+i < k pour 
i = 1, 2H, on H est le plus grand entier dans A/u. Les determinants sont calcules 
a partir des matrices (if + 1) x (if + 1) d'elements {h}-j et H x H d'elements {v}^j 
dont : 



Q*{k, AA, cj/A) = 1 - Prob{S{Lj) > k) 



(4.3) 



2i-l 





2j-2 





s=2i 



2i 





s=2j 



2i-i 

rus + {j - i)k I < i < j < H 

s=2i+l 
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En utilisant les definitions de V, hij et fy, nous avons pour k > 2 et u < A : 

Q*{k,\A,uj/A) = ^R*det\l/h,j\\det\l/vij\\ (4.4) 

V 

Dans la formule 14.41 : 

R* = m d^' (^ - df~^'p{N, AA) (4.5) 

H 

M = J2m2j+i (4.6) 

i=o 

d = 1 — uH/A et p(A^, AA) est la probabilite de poisson de N evenements et de 
rapport moyen AA. 

4.1.2 Approximation de Naus 

Une approximation tres utile de I'equation (14.41) a ete introduite par Naus 
[66] [65], basee sur les valeurs exactes des probabilites : 



Q2 = Q%k, 2^,1/2) 

et 

Qs = Q*ik, 3z^,l/3) 

on ip = Xuj et L = A/uo. 
On obtient alors [61] : 

Q*{k,,ljLML) = Ql[QllQt-^ (4.7) 



ou 



Ql = [F{k-l,i;)f -{k-l)p{k,i;)p{k-2,ij) 

-{k-l-ip)p{k,i;)F{k-3,ip) (4.8) 
Ql = [F{k-l,^)f-A^ + A2 + A,-A, (4.9) 
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et : 



A2 



A 



1 



2 p{k, i))F{k - 1, V) {(fc - 1) F{k -2,iP)- iPF{k -3,ip)} 
0.5 [p{k, {{k-l){k- 2) F{k - 3, ip) 
-2 {k - 2) '^F{k -4:,ip)+ ij^F{k - 5, V)} 



fe-i 



^3 



5^p(2A;-r,V)[F(r-l,V')]' 



fc-1 



^4 



^p(2A; - r, V)p(r, ^) {(r - 1) F(r - 2, ^Z^) - 7/'F(r - 3, ^Z-)} 



r=2 



Dans les formules au dessus F{k, ip) represente la distribution cumulative : 



et F{k, ?/') = pour < 0. 

Dans ce qui suit, nous utiliserons la methode scan statistics basee sur I'approxi- 
mation (14. 7p . 

4.1.3 Methode d'analyse 

Pour chercher I'existence d'une composante qui emet des rayons cosmiques dans 
une direction privilegiee, nous proposons d'appliquer I'approche scan statistics sur 
le flux de muons detectes par MACRO. Cette methode nous permet de chercher 
I'apparition de groupements d'evenements qui peuvent etre candidats de sources 
emettrices des rayons cosmiques dans I'espace environnent [67] . 
Les donnees utilisees pour cette analyse ont ete collectees a I'aide des tubes a strea- 
mer des 6 supermodules des parties inferieure et superieure du detecteur MACRO. 
Ce qui nous a permis de selectionner les muons reconstruits avec les 14 plans des 
tubes a streamer. 




fc 



(4.10) 



i=0 
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Donnees et criteres de selection 

Vu la nature de cette analyse, nous avons adopte des criteres de selection moins 
rigoureux par rapport a I'analyse de la distribution des temps d'arrivee des muons : 

- La duree des runs doit etre superieure a 2 heures (la valeur moyenne de la 
duree des runs est d'environ 5,5 heures). 

- La frequence des muons R est choisie dans I'intervalle 880 < R < 1006 fi/h 
et situee dans un intervalle ±2cr autour de la valeur moyenne. 

- Un temps mort d'acquisition inferieure a 2.5%. 

- L'efficacite de detection des tubes a streamer doit etre superieure a 90% et 
70% pour les fils et les strips respectivement. 

Un lot de 6411 runs a ete selectionne contenant un total d'evenements d'environ 
3.29 10^ muons. 

Methode d'analyse : 

Soit [Ai, Bi] I'intervalle de temps separant le temps du debut et de fin du run 
pour chaque run i. Une fenetre de temps de largeur uj balaye I'intervalle [AjjiJj] 
(voir figure \4A\\ . La Distribution du nombre moyen d'evenements dans une fenetre 
de largeur u dans chaque run est represente sur la figure 04. 2p pour u = 15 min, 
UJ = 5 min, uj = 1 min et = 30 s. 

Soit ki le nombre maximum d'evenements comptes durant le balayage. Sur les fi- 
gures (14. 3p et (14.41) est represente le nombre maximum d'evenements comptes durant 
le balayage des runs S{uj). 

Finalement, pour chaque run, on a calcule la probabilite pour que les fiuctuations 
statistiques puissent produire un cluster d'evenements aussi large que k^. Dans cette 
analyse le parametre uj est libre ; ici il est fixe au valeurs (15 min, 5 min, 1 min et 
30 s). Pour chaque run analyse, la distribution des probabilites scan statistics est 
representee sur les figures (14.51) et (14. 6p pour les differentes valeurs des fenetres de 
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Feiietie de lar^eur 

Fig. 4.1 Shematisation du balayage d'un run par une fenetre de largeur 
w 

temps {u = 30s , u = Imin) et {u = 5min , u = 15mm) respectivement. 

Pour les fenetres {w = Imin, 5mm, et 15mm), les figures (14.51) et (14. 6p montrent 
la presence d'un groupement d'evenements avec une faible valeur de probabilite scan 
statistics (inferieure a 10"'') identifies sur les runs 11056 et 11079. La position de ces 
groupements d'evenements est reperee a la fin de ces runs, ce qui ne permet pas de 
tirer une conclusion definitive sur la nature de I'exces. L'inspection de I'allure de la 
variation du nombre d'evenements en fonction du temps d'arrivee ne montre aucune 
divergence apparente ( voir figure (14.71) ). Les fieches sur la figure (14. 7p indiquent la 
position des groupements d'evenements. 
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Fig. 4.2 Distribution du nombre moyen d'evenements dans une fenetre 

de largeur ui dans chaque run pour o; = 15 min, uj = 5 min, ui = 1 min et 
a; = 30 s. 
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Fig. 4.3 Distributions du nombre maximum d'evenements S{u!) dans les 
fenetres de temps w = 30 s et w = 1 min. 
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Fig. 4.4 Distributions du nombre maximum d'evenements S{u!) dans 
les fenetres de temps w = 5 min et w = 15 min. On constate la pre- 
sence des fenetres avec un exces d'evenements (k=162 pour w= 5 min 
et k=346 pour w= 15 min) par rapport a la moyenne (<k>5min=i07 et 
<k>i5^in=277). 
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Fig. 4.5 Distributions de la probabilite scan statistics des runs selec- 
tionnes pour les valeurs des fenetres de temps w = 30 s et w = 1 min. 
Pour w=lmin, on note I'apparition d'une fenetre avec una faible valeur 
de probabilite scan statistics inferieure a 10~^. Cette fenetre correspond 
au run 11056. 
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Fig. 4.6 Distributions de la probabilite scan statistics des runs selec- 
tionnes pour les valeurs des fenetres de temps w = 5 min et w = 15 min. 
Pour w= 5 min et w= 15 min, on note I'apparition des fenetres avec 
de faibles valeurs de probabilite scan statistics inferieure a 10""^. Ces 
fenetres correspondent aux run 11056 et 11079. 
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Fig. 4.7 Xe nombre d'evenements en fonction du temps d'arrivee des 
muons pour les runs 11056 (a) et 11079 (b). Les fleches indiquent la 
position des groupements d'evenements reperes dans ces runs. 
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4.1.4 Conclusion 



L'analyse scan statistics portee sur 6411 runs selectionnees ne presente aucune 
deviation significative de I'hypothese nulle. Seules 2 runs ont ete identifies avec 
une faible probabilite dont I'exces a ete repere vers la fin du run, ceci pent etre 
attribue aux problemes d'acquisition des donnees qui peuvent parfois apparaitre 
vers la fermeture des runs. Ceci nous permet de confirmer I'absence d'eventuel 
cluster d'evenements dans le flux de muons et done aucune variation du flux n'a 
ete identifiee. 

4.2 Recherche de variations periodiques dans le flux 
de muons de MACRO 

4.2.1 Introduction 

Afin de chercher I'existence de variations periodiques dans le flux de muons de- 
tecte par MACRO, qui peuvent etre la signature d'une composante emettrice de 
rayons cosmiques. On se propose, dans ce paragraphe, d'analyser les signaux perio- 
diques dans les series de temps d'arrivee des muons. 

La methode de la transformee de Fourrier rapide (FFT) est la plus adaptee pour ce 
genre d'etude. Cette technique represente I'une des plus puissantes methodes qui 
permet d'identifier les fluctuations regulieres dans les series de temps. Cependant, 
une telle technique s'applique a des donnees echantionnees sur des intervalles uni- 
formement distribues et ne tient pas compte des bins vides. 

Dans le cas d'un echantillon de donnees non reguliers ou les sauts de donnees ap- 
paraissent, la FFT est appliquee apres un peuplement des sauts par les techniques 
d'interpolation, ce qui dissimule generalement les informations originales du signal. 
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Les methodes d'interpolations ont ete etudiees dans [70], il a ete recommande I'uti- 
lisation de la methode spectrale Lomb-Scargle [71] [72] qui permet de palier I'effet 
introduit par la presence de bins vides. 



4.2.2 Methode de Lomb-Scargle 

Les pannes on les periodes de maintenance du detecteur induisent des interrup- 
tions irregulieres dans la prise des donnees. La methode spectrale de Lomb-Scargle 
[7T][72| a ete developpee pour surpasser cet effet et conduire ainsi a la determina- 
tion de periodicites temporelles dans le flux de muons. 

Etant donnee N points mesures, la moyenne h et la variance a des donnees sont 
donnees par : 

N N 

h^-Y.^. <y'^j^^Y.^h^-hf (4.11) 

i=l j=l 

Le periodogramme de Lomb normalise est donne par : 



PnM = TT^ < ^ + 



2 



1 2 



2(T^ I cos^ ioitj — r) si'n? uj{tj — r) 

(4.12) 

oj = 2Txf > et r est defini par : 

y] • sin(2ujtj) 

tan{2ujT) = ^ (4.13) 

La constante r est un type d'excentrage qui rend Pn{w) invariant par translation de 
tons les ti par une constante. Lomb montre que ce type d'excentrage rend I'equation 
( 14.12P identique a I'equation que nous obtiendrons en estimant un ensemble de 
donnees harmoniques, a une frequence donnee tu, par un flt de moindre carree au 
modele : 

h{t) = A cosut + B sinut (4.14) 
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La methode de Lomb exploite les informations de chaque point contrairement a la 
FFT qui traite les intervalles. Un point tres frequent est que les points mesures hi 
sont la somme d'un signal periodique et le bruit gaussien independant. 
La presence d'un pic dans le spectre de P/v(co') pent reveler la presence d'une modu- 
lation an niveau du signal etudie. Dans notre cas, I'hypothese 'non nulle' considere 
que les donnees sont les valeurs independantes d'une gaussienne. 
Scargle montre [72] qu'a chaque uj particulier et dans le cas de Thypothese 'non nul- 
le', Pn{^) a une distribution de probabilite exponentielle avec une moyenne egale 
a 1. Autrement dit, la probabilite que P/v(co') sera entre les positifs z et z + dz est 
exp{—z)dz. 

II suit que, si nous balayons quelques M frequences independantes, la probabilite 
qu'il ne donne pas de valeurs plus grandes que z est (1 — e^)*^, ainsi : 

p{> z) = l-{l- e-' f (4.15) 

est la 'fausse alarme' de I'hypothese 'non nulle', qui est, le niveau de signification de 
n'importe quel pic observe dans Pn{w). Une petite valeur de la probabilite 'fausse 
alarme' indique une grande signification d'un signal periodique. 

4.2.3 Recherche de modulations periodiques dans le flux de 
muons 

En I'absence de modulations dans les mecanismes de production et de propaga- 
tion, le flux de rayons cosmiques, est distribue aleatoirement dans le temps. 
L'etude des variations journaliere du temps solaire et sideral du flux des particules 
montre que le detecteur presente une haute sensibilite pour l'etude des tres petites 
variations du flux de muons (de I'ordre de 10~^)[55]. Elle montre bien la grande 
sensibilite de MACRO meme aux faibles variations du flux de muons. Sur la figure 
(I3.4[ b) on observe bien la presence d'une variation periodique due a la variation 
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saisonniere du flux [54] . 

Notre objectif est d'analyser par la methode Lomb-Scargle les donnees de muons 
afin de chercher les modulations periodiques qui peuvent apparaitre. 
Nous avons echantillonne les muons sur des intervalles de temps de 15 min. Les 
bins deviant par plus de 3a par rapport au taux moyen mensuel ont ete elimines. 
Le nombre total des bins utilise est 160242 correspondant a 58% de la totalite de 
notre echantillon. 

Les resultats de I'analyse Lomb-Scargle appliquee sur notre echantillon de donnees 
sont representes sur la flgure (14.81) . 

Nous comparons le spectre obtenu par les donnees experimentales a celui de la si- 
mulation Monte Carlo. Dans cette derniere, nous avons ajouter des perturbations 
de meme niveau que les donnees reelles ( a savoir les ondes saisonniere, journaliere 
et siderale) et des intervalles de temps choisit ont ete distribues selon la sequence 
des series originales [54] . 

Sur la flgure ( 14.81) on remarque le grand pic a ~ 0.0027, il correspond a la variation 
saisonniere du flux. La flgure (14. 9p montre la region de frequence autour de la fre- 
quence journaliere solaire ou nous avons indique les frequences qui correspondent 
aux ondes siderale et anti-siderale, pour les donnees reelles (flgure en haut) et pour 
la simulation Monte Carlo (flgure en bas). Afln d'eliminer le bruit dans le spectre 
des puissances, nous avons considere un echantillonnage, aussi bien pour les don- 
nees reelles que pour les donnees simulees, selon la formule suivante : 

ou Ni est le contenu des bins originales et iV(Ar) est le contenu moyen dans I'in- 
tervalle de temps ±Ar (Ar a ete choisit egale 1 jour). Le pic a la frequence jour~^ 
a une signiflcation statistique de I'ordre de ~ 2.3 a. Le signal correspondant a la 
variation siderale est observe, mais des pics avec une taille similaire (ou encore plus 
large) sont aussi presents dans le spectre. 
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Fig. 4.8 En haut, la puissance spectrale de Lomb en fonction de la fre- 
quence [jour^] pour les donnees experimentales. Le pic a ~ 0.0027 cor- 
respond a la variation saisonniere. En has, nous presentons les resultats 
de la simulation en considerant un echantillon qui contient des periodi- 
cites de mimes ordre de grandeur que les donnees reelles (saisonniere, 
journaliere et siderale)[5^. 
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Conclure que les ondes journalieres solaire et siderale sont reelles depend essentiel- 
lement de la stabilite de leur amplitudes et leur phases avec le temps. 
Apres I'analyse des series de temps avec la methode de Lomb-Scargle, nous avons 
trouve a nouveau la modulation saisonniere (figure 14.81) . et des signaux aux po- 
sitions des ondes journalieres solaire et siderale (figure US]) . Aucune deviation de 
I'hypothese nulle n'a ete signale. 

4.3 Conclusion 

Pour I'analyse des series de temps des muons, detectes par MACRO, deux ap- 
proches complementaires ont ete considerees : la recherche de periodicite en utilisant 
la methode Lomb-Scargle et la recherche de cluster d'evenement en utilisant la me- 
thode scan statistics. Les deux methodes completent I'analyse de la distribution 
des temps d'arrivee des muons de haute energie menee dans le chapitre[3l Pour les 
deux approches aucune deviation de "I'hypothese nulle" n'a ete mise en evidence. 
Ce resultat confirme les conclusions du chapitre [3] et permet a nouveau de conclure 
sur la nature aleatoire des temps d'arrivee des muons. 
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Fig. 4.9 La region de frequence autour de Vonde journaliere pour les 
donnees reelles en (haut) et pour la simulation Monte Carlo (en has). 
Les fleches montrent la position des pics journalier solaire, sideral et 
anti-sederal. 
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Chapitre 5 



Etude de la perte d'energie des 
nuclearites 



5.1 introduction 

Proposee par Witten [73], la matiere nucleaire etrange "Strange Quark Matter" 
(SQM) presente un etat de la chromodynamique quantique. C'est un agregat de 
quarks u, d et s en nombre egale, entoure par un nuage d'electrons assurant la 
stabilite electrique et formant ainsi la dite "nuclearite". 

Un grand interet a ete apporte a I'etude des nuclearites, Bodmer [74] a suggere 
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que les nuclearites peuvent etre plus stables que la matiere nucleaire ordinaire. Farhi 
et Jaffe [75] out montre que les nuclearites peuvent etre stables sur un large domaine 
de masses, allant de celles des noyaux les plus legers (quelques quarks) jusqu'a celles 
des etoiles a neutrons (~ 10^^ quarks), a conditions qu'ils se composent de quarks 
u, d et s en quantites approximativement egales. En effet, la taille des nuclearites 
pent depasser largement celle des noyaux les plus massifs, car contrairement a ces 
derniers, il n'y pas de limitation due a la repulsion coulombienne. 

En partant de ces conditions, la matiere etrange doit etre extremement stable. 
Elle pent former des systemes lies avec n'importe quel nombre de quarks, et pent 
etre consideree comme un candidat a la matiere obscure de I'univers [76] . II est 
suppose que les nuclearites out ete produits aux premiers instants de I'univers, 
(T > 200MeV ett< lO^s) [77]. La production des nuclearites, pent avoir lieu aussi 
dans les phenomenes astrophysiques violents se produisant au sein de I'univers dont 
on pent citer les collisions entre etoiles a neutrons et les explosions des etoiles [78] . 

La possibilite que les nuclearites bombardent la terre a ete introduite par de 
Rujula et Glashow [79] . 

Dans ce qui suit, nous discuterons la perte d'energie de ces particules et leur 
possibilite pour arriver aux differents niveaux de detection. 

5.2 Proprietes des nuclearites 

On pent approximer un nuclearite par une sphere de rayon R, dans laquelle les 
quarks u, d et s sont confines dans un volume V avec une pression B dite pression 
de "bag". Les nuclearites sont decrits par le "bag model" (Modele du Sac) [77] 
dans lequel les quarks representent un gaz de Fermi degenere. La matiere nucleaire 
(MN), pent etre vue comme une association de plusieurs sacs maintenus ensemble 
par une energie de liaison (a titre d'exemple le deuton se compose d'un sac uud 
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et d'un sac udd). Cette idee peut etre etendue a des sacs contenant plus de trois 
quarks dans un meme sac, c'est la matiere de quarks |73| . Lorsque la matiere de 
quarks contient des quarks "strange" c'est la matiere nucleaire etrange "SQM". 
II est suppose que la charge est distribuee de maniere uniforme a I'interieur du 
nuclearite et que I'equilibre "chimique" entre les quarks u, d, s et les electrons, est 
maintenu par les interactions faibles |77j 



d ^ u + e + Pe 

u + e —>■ d + 
s — >• M + e + z7g 

u + e ^ s + Ue 
s + u < — »• u + d 

L'etat du nuclearite est determine par le potentiel thermodynamique = u, d, s). 
Ce dernier est fonction du potentiel chimique /ij, de la masse du quark s(ms), de 
la pression de bag B et du volume V. 

Les densites des especes sont donnees par : 



et la densite baryonique 



n. = -f^ (6.1) 



n=— ^ (5.2) 



Afin d'assurer I'equilibre du potentiel chimique des especes de quarks, SQM doit 
avoir un nombre de quarks s un peu moins que le nombre des quarks u et d [77] . 
Cependant le noyau du nuclearite doit avoir une charge electrique positive com- 
pensee par le nombre d'electron du cortege electronique = {Nd + A^s)/3, ou N^, 
Ng et Ne sont les nombres de quarks d, s et les electrons respectivement ( il a ete 



^On appelle noyau du nuclearite le systeme quarks sans cortege electronique 
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considere que Nd = Nu) [80] . 

Soit Rn le rayon du noyau du nuclearite, le systeme nuclearite aura un rayon 
constant ~ 1 A pour 

Pour Rn > 10^ fm, tons les electrons sont a I'interieure du sac de quark. 

Pour 10^ < Rn < 10^ fm, une fraction des electrons est a I'interieure du sac et une 

autre a I'exterieure, il donne an nuclearite dans ce cas une dimension globale de 

~ 10~^ cm = 10^ fm. Dans ces conditions le nuclearite est similaire a I'atome de 

Bohr. 

La figure (15.11) illustre une distribution spatiale du systeme noyau+electrons. Notons 
que pour Rn > 10^ fm, le nuclearite est considere comme un atome de Thomson 

m- 




Rj^=10^fm R^^lO^fm RN=10^fm Rj^^lO^fm Rn= 10^ fm 

Fig. 5.1 Dimensions du sac de quarks (Rn) et le systeme 
noyau+electrons (nuclearite). Les points indiquent les electrons, le bord 
du sac de quarks est indique par la ligne solide ; le bord du systeme 
noyau+ electron pour les faibles masses est indique par la ligne en poin- 
tilles. 
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Pour les nuclearites de masses superieures a 1.5-10 g — 10 GeV/c , la relation 
entre la masse et le rayon prend la forme suivante (masse oc volume, Mat oc Vat oc 



R 



■N 



3M, 



N 



AnpN 



1/3 



(5.3) 



Pour Mm = 1.5 -10 ^ g, le rayon du nuclearite devient : 



R 



TV 



1/3 



IQ-^cm = 1 A 



(5.4) 



^47r3.5 ■ 1014^ 

Pour Mm = 10^^ GeV/c^, R^ - 10^ fm (1000)^/^ ^ iqg 

Le rayon du noyau du nuclearite -Rat en fonction de sa masse est represents 
sur la figure (15. 2p . 
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Fig. 5.2 Rayon du nuclearite en fonction de la masse (modele du sac de 
quarks) : ligne Solide. La ligne en pointilles donne le rayon du systeme 
noyau+ electrons (nuclearite). Nous indiquons aussi les regions d'acces- 
sibilite en terme de masse par rapport aux experiences MACRO et SLIM. 
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5.3 Perte d'energie des nuclearites 



Dans la plupart des travaux qui calculent la perte d'energie des nuclearites dans 
divers milieux detecteurs, on ne tient pas compte de la contribution de I'atmosphere. 
Cette derniere pent etre considerable dans le cas des experiences souterraines ou 
celles situees au niveau de la mer vu que la densite de I'atmosphere augmente tout 
en se rapprochant de la mer. 

Dans notre travail, on a introduit la contribution de la perte d'energie dans 
I'atmosphere afin de determiner la vitesse minimale avec laquelle les nuclearites 
peuvent atteindre differents niveaux ou sont installes les detecteurs, ainsi que la 
vitesse d'arrivee des nuclearites en considerant qu'ils atteignent la haute atmosphere 
avec une vitesse P = 10~^. 

La perte d'energie des nuclearites traversant la matiere est due principalement 
aux collisions elastiques ou quasi-elastiques avec les atomes du milieu [79j. Aux 
faibles vitesses, I'interaction des nuclearites se fait par des collisions elastiques avec 
les molecules du milieu. 

Une partie des objets galactiques massifs non lumineux pent exister sous forme 
de nuclearites [79] . 

La densite locale de la matiere obscure est de I'ordre de 10'"^^ g/cm^ avec des 
vitesses d'echappement de la galaxie avoisinant 250 km/s. Le rapport de perte 
d'energie des nuclearites ayant une telle vitesse est p9] : 



Ou A est la section efficace effective du nuclearite, v sa vitesse et p la densite du 
milieu traverse. 

La section efficace geometrique des nuclearites est donnee par : 



dE , . 
— = -Apv^ 
ax 



(5.5) 



A = 7iR'^ 



(5.6) 
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Comme il a ete discute au dessus, I'aire effective des nuclearites de faibles dimensions 
est controlee par son nuage electronique qui ne pent jamais etre inferieur a ~ 1 A. 
A faibles masses {A < 10''), le nuage electronique n'est pas a I'interieur du nuclearite 
[75] . Pour des masses (10^ < ^4 < 10^^) une partie du nuage electronique va se 
trouver a I'interieur. 

Pour des masses A > 10^^, tout le nuage electronique sera contenu a I'interieur du 
nuclearite, ce dernier devient neutre (le nuage electronique fait ecran aux charges 
positives internes de la particule). 
La section effective d'interaction est donnee par [79] : 



„ , TT X 10 pour M < l.bnq 

A cm" <^ 5.7 
7r(^)2/3 pour M > l.5ng 

Ou Pn est la densite de la matiere nucleaire etrange, estimee etre egale a ~ 3.6 x 
10^^ g/cm^ [81]. 

La section efficace d'interaction des nuclearites en fonction de leur masse est 
representee sur la figure (15.31) . 

Selon I'equation (15.51) . la vitesse des nuclearites diminue de fagon exponentielle, 
lors de la traversee de la matiere selon I'equation suivante : 



v{L) = vq exp 



A 
M 



L 



p dx 



(5.8) 



Oil L est la longueur traversee et vo la vitesse d'entree dans I'atmosphere terrestre. 
Ainsi le parcours R des nuclearites pent s'ecrire sous la forme : 

Oil Vc est la vitesse minimale pour traverser la longueur L. 
Pour une vitesse P ~ 10"^, le parcours est donnee par : 

3 X 10MM/ln(7]i/3 pour M> 1.5ng 

R[g/cm ) = < (5.10) 
2.3 X 10^ [M/lng]^/^ pour M < l.bng 
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masse des nucleorites (g) 



Fig. 5.3 Section efficace d'interaction des nuclearites en fonction de 
leur masse. 
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Masse des nuclearites (g) 



Fig. 5.4 Parcours des nuclearites en fonction de leur masse. 

Sur la figure (15.41) on presente le parcours en fonction de la masse des nuclearites. 
II est a noter que les nuclearites ayant des masses M > 4.5 x lO^^^g traversent 
I'atmosphere et arrivent a la surface de la terre, alors que ceux ayant des masses 
M >0.1 g peuvent traverser la terre sans etre arretes. 

5.4 Perte d'energie des nuclearites dans la terre 

5.4.1 Le modele interne de la terre 

La composition exacte de la matiere qui constitue la terre presente de grandes 
controverses [37], par consequent la perte d'energie dans la terre ne pent etre qu'ap- 
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Fig. 5.5 Profil de la densite terrestre est represente par la courbe discon- 
tinue (voir Ref. [37^); La courbe en pointilles represente une approxi- 
mation de la densite terrestre (tq = rayon de la terre). 

proximative. 

Le manteau terrestre devrait etre forme de dioxyde de silicium Si02. Le profil 
de la densite terrestre, que nous avons adopte pour notre etude, est represente 
en pointilles sur la figure (15.51) . ou Ton pent observer trois parties : le noyau, le 
manteau et la croiite. Un modele simple a ete applique, dans lequel la densite et la 
composition de chaque partie est uniforme et elle est presentee en lignes continues 
sur la fig. (15. 5p [82] . Dans ce modele le noyau est constitue de fer, avec une densite 
Pfer = ll.5g/cm^ et une conductivite de 1.6 x lO^^s"^ ; le manteau est constitue de 
silicium avec une densite psi = i.Sg/cm. Le rayon du noyau constitue 0.54 fois le 
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rayon terrestre. 

Les nuclearites perdent leur energie dans la terre par des collisions elastiques avec 
les atomes du milieu. On represente la perte d'energie des nulearites de differentes 
masses dans le noyau de la terre sur les figures (I5.6p et dans le manteau terrestre 




S -4 -3 

10 10 10 

^=v/c 



Fig. 5.6 Perte d'energie des nuclearites de differentes masses dans le 
noyau de la terre. 



5.5 Perte d'energie des nuclearites dans I'atmosphere 

La complexite du calcul de la perte d'energie des nulcearites dans I'atmosphere 
vient de la variation de la densite de I'atmosphere avec la profondeur. Ainsi la 
perte d'energie depend de cette derniere. Le calcul de la perte d'energie se fait 
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Fig. 5.7 Perte d'energie des nuclearites de differentes masses dans le 
manteau terrestre. 

en se basant sur le fait que I'atmosphere est divisee en succession de couches de 
differentes epaisseurs et de densites. 

Dans ce travail la densite de Tatmosphere a ete calculee en utilisant la parame- 
trisation de Shibata [83] : 



ou les constantes a = 1.2 x l{)~^g.cm~^ et 6 ~ 8.57 x lO^cm, H est I'altitude totale 
de Tatmosphere (~ 50 km) et L la profondeur de penetration dans I'atmosphere. 
La variation de la densite de I'atmosphere en fonction de I'altitude est representee 
sur la figure (15. 8p . 



La perte d'energie des nuclearites dans I'atmosphere est montree sur la figure (15.90 
pour les masses M = 8.4 x 10^^ GeV/c^ et M < 8.4 x 10^^ GeV/c^. 




p{h) = ae b = ae >> 



(5.11) 
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Fig. 5.8 Variation de la densite de V atmosphere en fonction de I'alti- 
tude. 




profondeur dans I atmosphere (km) 



Fig. 5.9 Perte d'energie dans I'atmosphere des nuclearites de masse 
M < 8.4 lO^^GeV/c^ (a) et M ^ 8.4 lO^^GeF/c^ (b). 
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5.5.1 La region d'accessibilite des nuclearites 



Pour calculer la region d'accessibilite dans le plan (masse, /?), I'integration de 
I'equation ( 15.5p a ete resolue analytiquement : 



A partir de I'equation ( IS.Sp on determine la vitesse minimale que pent avoir les 
nuclearites de differentes masses pour atteindre les differents niveaux dans I'atmo- 
sphere. 

Sur la figue (15.101) on montre la vitesse des nuclearites avec laquelle les differentes 
masses peuvent atteindre les niveaux correspondants a : 

- I'experience de ballon CAKE (40 km) ^7], 

- au possibles experiences portees sur des avions civils (11km), 

- au laboratoire de Chacaltaya (I'experience SLIM, 5.29 km) [84] [85]. 

- au niveau de la mer, 

- I'experience MACRO [86] est aussi indue. Le seuil de detection dans le CR39 
(correspondant a une perte d'energie restreinte REL= 200 MeV cm^), et 
dans le MACROFOL (REL=2500 MeV g'^ cm^) est represente. 

Nous avons considere les nuclearites ayant une vitesse initial (3 = 10~^ en haut de 
I'atmosphere. Un calcul pareil a ete fait pour /3 = 10~^, 10""^, 10^^ et n'a montre 
aucun effet de la vitesse initial sur le resultat obtenu. 

La diminution de la vitesse seuil pour les nuclearites ayant une masse superieure a 
8.4 X lO^'* GeV est due au changement de la section efficace selon I'equation (15. 7p . 
Dans une experience portee a I'altitude de Chacaltaya la masse minimale des nu- 
clearites detectables s'abaisse d'un facteur 2 par rapport a d'autres experiences 
installees au niveau de la mer. 




(5.12) 
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Fig. 5.10 Vitesse d'arrivee des nuclearites a differentes profondeurs 
en fonction de la masse, en considerant que la vitesse d'arrivee a la 
haute atmosphere (3 = 10 : les lignes solides. Les lignes en pointilles 
representent le seuil de detection dans le CR39 et dans le MAKROFOL. 

Les conditions de detection des nuclearites dans le CR39 sont exprimees par 
la vitesse d'entree minimale en haut de I'atmosphere (50 Km) en fonction de la 
masse des nuclearites. Pour les differentes locations d'experiences, on represente 
sur la figure (15. lip la vitesse d'entree minimale en fonction de la masse. Dans ce 
cas, la contrainte est que le nuclearite possede une vitesse minimale au niveau du 
detecteur afin qu'il puisse produire une trace ; notons que pour toutes les experiences 
le seuil REL de detection est pris egale a 200 MeV cm? pour le CR39 et 
250 MeV g~^ cm? pour le MAKROFOL. Le changement de I'allure de la courbe 
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est du a la variation de la section efficace en fonction de la masse selon I'equation 




Masse des nuclearites (Gev/c^) 



Fig. 5.11 Conditions de detection dans le CR39 et le MAKROFOL pour 
differentes experiences situees a diverses altitudes. 



122 



5.6 Conclusion 



En se basant sur la formule de De Rujula (15.51) relative a la perte d'energie des 
nuclearites et en utilisant la modelisation de Sliibata pour le calcul de la densite 
de ratmosphere, nous avons reussi a evaluer la perte d'energie des nuclearites le 
long de leur parcours dans Tatmosphere. Pour les differentes experiences realisees 
a differentes altitudes dans Tatmosphere (CAKE (40km), experiences portees sur 
des avions civil (11km), SLIM (5.29km), experiences situees au niveau de la mer) 
nous avons determine la perte d'energie dans le plan des parametres (masse,/?). En 
comparant le resultat obtenu dans la figure (I5.10p avec celui du [87], on remarque 
une difference de masse d'un facteur de ~ 3. Ceci est du a la fagon avec laquelle a 
ete calculee la densite de I'atmosphere. 

Les nuclearites ayant une masse > 7 lO^^GeV/c^ ne sentent pas la presence de I'at- 
mosphere. Au dessous de cette masse, et pour une experience portee a une altitude 
de 5.29 Km (I'experience SLIM ), le minimum de masse detectable decroit d'un 
facteur 2 par rapport aux experiences installees au niveau de la mer. 
Pour I'experience souterraine MACRO (3800 m.w.e), les nuclearites de masse su- 
perieure a 1.53 lO^^GeV/c^ (figure ISTTOll ne subissent aucune perte d'energie, leur 
Vitesse /3 est intacte. Par contre le parametre /3 decroit exponentiellement pour des 
masses dans I'intervalle [1.6 10^^, Q 10^^]GeV/c^. Le CR39 du detecteur MACRO 
est sensible aux nuclearites de masses > 2 lO^^GeV/c^ pour lesquelles 3.37 10^"^ < 
(3 < 0.9 10"^. Dans le cas du MAKROFOL, MACRO est sensible aux nuclearites 
de masses > 3 lO^^GeV/c^ pour lesquelles 0.12 10"^ < f3 < 0.9 10'^) . 
La contribution de I'atmosphere dans le calcul de la perte d'energie devient consi- 
derable pour les nuclearites avec des masses < 7 10^'^GeV/c^ et ayant une vitesse 
P = 10^^ en haut de I'atmosphere. 
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Conclusion 



Le travail que nous avons mene au pres de I'experience MACRO s'est articule 
autour de la recherche des correlations possibles dans le temps d'arrivee des muons 
cosmiques ainsi que sur la recherche des variations dans le flux de muons. 

Nous disposons d'un lot de donnees d'environ 38 millions de muons d'energie 
superieure a 1.3 TeV au sommet de la montagne du Gran Sasso, collecte a I'aide 
des tubes a streamer du detecteur MACRO dans sa configuration complete. 

L'analyse des distributions des temps d'arrivees des muons singuliers, doubles 
et multiples arrivant de toutes les directions du ciel ainsi que ceux venant des zones 
selectionnes formant des cones, montre qu'elles sont compatibles avec la fonction 
gamma. Ce qui confirme I'hypothese d'arrivee aleatoire des RC. Ce resultat a ete 
appuye par le test de Kolmogorov-Smirnov, dont aucune deviation entre la distri- 
bution theorique et celle experimentale n'a ete mise en evidence. 
Nos resultats ne sont pas ainsi analogue a ceux obtenus par Bath et al. [52], ou par 
Badino et al. [53j, et ceci pour des energies de I'ordre de 20 TeV (// singuliers) et 
superieures a 20TeV {fi doubles et multiples). 

La recherche de cluster d'evenements, dans le flux de muons, par la methode 
scan statistics n'a presente aucune deviation significative de I'hypothese nuUe. Ce 
qui nous a permis de confirmer I'absence d'eventuel cluster d'evenements dans le 
fiux de muons et done aucune variation du fiux n'a ete signale. 



124 



En analysant les series de temps des muons par la recherche des periodicites avec la 
methode Lomb-Scargle, nous avons observe la modulation saisonniere et des signaux 
a la position des ondes journaliere solaire et siderale. Aucune autres deviation n'a 
ete signale. En se basant sur la formule de De Rujula (15 .Sp relative au calcul de 
la perte d'energie des nuclearites et en utilisant la modelisation de Shibata pour 
le calcul de la densite de I'atmosphere, la contribution de la perte d'energie dans 
I'atmosphere a ete introduite afin de determiner la vitesse minimale avec laquelle 
les nuclearites peuvent atteindre differents niveaux ou sont installes les detecteurs. 
La contribution de I'atmosphere, pour les nuclearites ayant une vitesse (3 = 10^^ 
en haut de I'atmosphere, devient considerable pour les masses < 7 lO^'^GeV/c^ . 
Au dessous de cette masse, et pour une experience portee a une altitude de 5.29 
Km (I'experience SLIM ), la masse minimale detectable decroit d'un facteur 2 
par rapport aux experiences installees au niveau de la mer. Pour les experiences 
souterraines, les masses superieures a 1.53 lO^^GeV/c^ ne subissent aucune perte 
d'energie. 



125 



Bibliographie 



[1] M. Ambrosio et al., Eur. Phys. J. C36 (2004) 323; 
M. Ambrosio et al., Phys. Lett. B566 (2003) 35. 

[2] M. Ambrosio et al., Phys. Lett. B478 (2000) 5. 

[3] M. Ambrosio et al., Phys. Lett. B434 (1998) 451 ; 
M. Ambrosio et al, Phys. Lett. B517 (2001) 59. 

[4] M. Ambrosio et al., Astrop. Phys. 19 (2003) 313; 

[5] D. Sivoukhine, cours de Physique generale. Editions Mir de Moscou. 

[6] K. Greisen, Phys. Rev. Lett. 16 (1966) 748 ; 

G.T. Zatsepin et V.A. Kuzmin, Prisma Zh. Exp. Theor. Phys. 4 (1966) 114. 

[7] J. R. Ellis et al. Phys. Lett. B 247 (1990) 257. 

[8] M. Takeda et al., Astrop. Phys. 19 (2003) 447. 

[9] R. U. Abasi et al., Phys. Rev. Lett. 92 (2004) 151101. 

[10] G. Matthiae, "Results from the Auger Observatory on high-energy cosmic 
rays", AIP Conf. Proc. 1059 (2008) 143. 

[11] P. Sokolsky, "Introduction to Ultrahigh Energy Cosmic Ray Physics", Addison- 
Wesley 1989. 

[12] J. A. Simpson, Ann. Rev. Nucl. Part. Sci. 33 (1983) 323. 

126 



[13] T. K Gaisser, "Cosmic Rays and Particle Physics", Cambridge University Press 
1990. 

[14] E. Fermi, Phys. Rev. 75 (1949) 1169. 

[15] M. S. Longair, "High Energy Astrophisics vol.1 et 2", Cambridge University 
Press 1992. 

[16] V. L. Ginzburg et S. I. Syrotskii, "The origin of Cosmic Rays", Oxford Perga- 
mon Press 1964. 

[17] M. A. Ruderman and P. G. Sutherland, Astrophys. J. 196 (1975) 51. 



[18] A. V. Olinto, "Ultra high energy cosmic ray accelerators", astro-ph/9911154 



[19] N. Hayashida et al., "The anisotropy of cosmic ray arrival direction arround 

lO^^eV", Astrop. Phys. Vol. 10 (1999) 303. 
[20] R. C. Hartman et al., Astrophys. J. 230 (1979) 579. 
[21] D. Bird et al., Ap. J. 511 (1999) 739. 
[22] M. Takeda et al., Ap. J. 522 (1999) 225. 
[23] M. Drees et al., Phys. Rev. D39 (1989) 1310. 

[24] H. Dekhissi, These de Doctorat d'Etat,"£^^MC?e spatio-temporelle des muons 

cosmiques avec le detecteur MACRO et etude du detecteur plastique CR39", 

Universite Mohamed Premier, Oujda (1994). 
[25] F. Maaroufi, These de Doctorat d^Etat" Etude de I'anisoptropie temporelle et 

spatiale des rayons cosmiques avec le detecteur MACRO", Universite Mohamed 

Premier, Oujda (1999). 
[26] S. Ahlen et al., Nucl. Phys. B370 (1992) 432. 

[27] R. Fantini, Tesi di Dottorato, "Ricerca di monopoli magnetici con i tubi a 

streamer di MACRO", Universita di Bologna (1996). 
[28] M. Sitta, Tesi di Dottorato, " Rivelazione di monopoli magnetici con i tubi 

a streamer dell' esperimento MACRO al Gran Sasso" , Universita di Torino 

(1994). 



127 



[29] S. Ahlen et al., Nucl. Instr. Meth. Phys. Res. A324 (1993) 337. 

[30] C. De Marzo et al., Nuovo Cimento C9 (1986) 281. 

[31] M.Ambrosio et al. Nucl. Instr. Meth. Phys. Res. A486 (2002) 663. 

[32] M. Ambrosio et al., Astrop. Phys. 6 (1997) 113. 

[33] M. Ambrosio et al., Astrop. Phys. 4 (1995) 33. 

[34] S. Cecchini et al., Nuovo Cimento 109 (1996) 1119. 

[35] S. Cecchini et al., Astrop. Phys. 1 (1993) 369. 

[36] M. Ambrosio et al., Phys. Rev. Lett. 72 (1994) 608. 

[37] D. Anderson, "Theory of the earth", Blackwell Scientific Publications (1989) ; 
J. P. Poirier, "Introduction to the physics of the earth's interior", Cambridge 
University Press (1991). 

[38] M. Ambrosio et al., Eur. Phys. J. C25 (2002) 511. 

[39] M. Ambrosio et al., Eur. Phys. J. C26 (2002) 163. 

[40] S. Ahlen et al., Phys. Rev. Lett. 72 (1994) 3046. 

[41] M. Ambrosio et al., Phys. Lett. B406 (1997) 249. 

[42] S. Ahlen et al., Astrop. Phys. 1 (1992) 11. 

[43] S. Ahlen et al., Astrophysical J. 412 (1993) 301. 

[44] N. Giglietto, Nucl. Phys. B, Proc. Suppl., Vol. 412 (1993) 301. 

[45] S. Ahlen et al., Phys. Rev. D46 (1992) 4836. 

[46] S. Ahlen et al., Phys. Rev. D50 (1994) 3046. 

[47] S. Ahlen et al., Phys. Rev. Lett. 69 (1992) 1860. 

[48] P. Ullio, "Dark Matter and Dark Energy", VII School on Non-Accelerator 
Astroparticle Physic (Trieste 2004). 

[49] T. C. Weekes, Nature 233 (1971) 129. 



128 



[50] C. Morello et al., Nuovo Cimento C7 (1985) 682. 

[51] P. Serri, Tesi di laurea, "studio dei muoni penetrante neU'esperimento NUSEX 
del Monte Bianco", Universita di Milano (1984). 

[52] C. L. Bhat et al., Nature 288 (1980) 146. 

[53] G. Badino et al., Nuovo Cimento 28 (1980) 93. 

[54] M. Ambrosio et al., Astrop. Phys. 7 (1997) 109. 

[55] M. Ambrosio et al., Phys. Rev. D67 (2003) 042002. 

[56] H. Press et al., "Numerical Recipes in C++", Cambridge University Press. 

[57] F. James, "Minuit Reference Manual" CERN program library entry D506, 
CERN Geneva (1994). 

[58] M. Ambrosio et al., Astrop. Phys. 18 (2003) 615. 

[59] W. T. Eadi et al., "Statistical methods in experimental physics" North-Holland 
(1971) 269. 

[60] Y. Becherini, H. Dekhissi, A. Moussa et al., Astrop. Phys. 23 (2005) 341. 

[61] F. Terranova, Nucl. Instr. Meth. Phys. Res. A519 (2004) 659. 

[62] J. Glaz, J. Naus and S. Wallenstein, "Scan Statistcis" Springer, New York 
2001. 

[63] K. J. Orford, J. Phys. G26 (2000) Rl. 
[64] R. J. Huntington and J. I. Naus, Ann. Prob. 3 (1975) 898. 
[65] S. Wallenstein and J. I. Naus, Ann. Prob. 1 (1973) 188. 
[66] J. I. Naus, J. Amer, Stat. Ass. 77 (1982) 177. 

[67] S. Cecchini, H. Dekhissi, A. Moussa et al., EPL 87 (2009) 39001, 
larXiv :0907.,5235l [hep-ex]. 

[68] S. E. Aim, "On the distribution of scan statistics in a Poisson process, in 
Probability and mathematical statistics" A. Gut and L. Heist eds., Upsalla 
University press 1983. 



129 



[69] K. Ackerstaffet al., Phys. Lett. B429 (1998) 399. 

[70] P. Laguna et al., "Power spectral density of unevenly sampled data by least- 
square analysis : performance and application to heart rate signals". IEEE 
Trans. Biomed. Eng. 1998. 

[71] N. R. Lomb, Astrophys. Space Sci. 39 (1976) 447. 

[72] J.D. Scargle, Ap. J. 263 (1982) 835. 

[73] E. Witten, Phys. Rev. D30 (1984) 272. 

[74] A. Bodmer, Pys. Rev. D4 (1971) 1601. 

[75] E. Farhi and R. L. Jaffe, Phys. Rev. D30 (1984) 2379. 

[76] J. H. Appelgate and C.J. Hogan, Phys. Rev. D31 (1985) 3037; 

A. Alcock and A. Olinto, Ann. Rev. Nucl. Part. Sci. 38 (1988) 161 ; 

J. Madsen, H. Heiselberg and K. Riisager, Phys. Rev. D34 (1986) 2947; 

J. Madsen and M.L. Olsen, Phys. Rev. D43 (1991) 1069. 

[77] M. Kasuya, T. Saito, M. Yasue Phys. Rev. D47 (1993) 2153. 

[78] P. Haensel, J.L. Zdunik and R. Schaeffer, Astron. Astrophys. 160 (1986) 121 ; 
A. Alcock, E. Farhi and A. Olinto, Astrophysical. J. 310 (1986) 261 ; 
D. Kaplan and A. Nelson, Phys. Lett. 157 (1986) 57; 

H.A. Bethe, G.E. Brown and J. Cooperstein, Nucl. Phys. A462 (1987) 791. 

[79] A. De Riijula, S.L. Glashow, Nature 312 (1984) 734 ; 
A. De Riijula, Nucl. Phys. A434 (1985) 605. 



[80] D. Bakari et al., hep-ex/0004019 



[81] E. Farhi and R. L. Jaffe, Phys. Rev. D32 (1985) 2452. 

[82] J. Hong, "Search for GUT magnetic monopoles and other supermassive par- 
ticles with The MACRO detector", California University 1993. 

[83] T. Shibata, Prog. Theor. Phys. 57 (1977) 882. 

[84] SLIM Collaboration, Nuovo Cim. Vol. 24 C, N. 4-5 April/May 2001. 



130 



[85] S. Cecchini et al., Eur. Phys. J. C57 (2008) 525; 
S. Balestra et al., Eur. Phys. J. C55 (2008) 57. 

[86] M. Ambrosio et al., Eur. Phys. J. C13 (2000) 453. 

[87] D. Bakari et al., |hep-ex/000 3028. 



131 



Resume 



La radiation cosmique est une composante principale de la galaxie, vu que sa 
densite d'energie est comparable, si elle n'est pas superieure, a celle des autres ra- 
diations presentes dans I'univers. Elle joue le role des messagers de mecanismes 
astropliysiques pouvant mettre en jeu des energies colossales. Elle est la seule ma- 
tiere qui arrive de I'exterieure de notre systeme solaire et elle constitue I'unique 
moyen pour I'etude des particules de tres haute energie allant jusqu'a 1019 eV. 

Parmi les questions de recherche qui presentent un interet fondamental, on dis- 
tingue I'etude de I'anisotropie temporelle du flux de muons cosmiques. Des mo- 
dulations dans les distributions temporelles du temps d'arrivee de muons ont ete 
observees par quelques experiences. Les resultats de ces experiences n'ont pas ete 
conflrmes par d'autres groupes et le probleme reste confus. 

A cet efl^et, la grande partie de cette these a ete consacree a I'etude du caractere 
aleatoire du temps d'arrivee des rayons cosmiques a travers I'etude des variations 
du flux des mouns cosmiques coUecte avec le detecteur souterrain MACRO. Deux 
techniques ont ete utilisees : 

- L'etude des distributions des temps d'arrivee des muons cosmiques. 

- L'etude des variations du flux de muons en utilisant deux approches : 

1. methode Scan-statics. 

2. methode de Lomb-Scargie. 

En s'appuyant sur les deux etudes cite au-dessus, la recherche de possible correla- 
tions temporelles on de periodicite dans le flux de muons cosmique de MACRO, 
nous a permis de conflrmer I'hypothese d'arrivee aleatoire des rayons cosmiques. 



La derniere partie de cette these a ete consacree a I'etude de la perte d'energie 
des nuclearites et leur possibilite pour arriver aux differents niveaux de detection 
dans I'atmosphere. La contribution de I'atmosphere dans la perte d'energie pour les 
nuclearites ayant une vitesse {3 = 10~^ en haut de I'atmosphere, devient considerable 
pour les masses < 7 Au dessous de cette masse, et pour une experience 

portee a une altitude de 5.29 Km (I'experience SLIM ), la masse minimale detec- 
table decroit d'un facteur 2 par rapport aux experiences installees au niveau de la 
mere. Pour les experiences souterraines, les masses superieures a 1.53 lO^^GeV/c^ 
ne subissent aucune perte d'energie. 
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